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The behavior of two individuals, consisting of effort which results 
in output, is considered to be determined by a satisfaction function which 
depends on remuneration (receiving part of the output) and on the effort 
expended... The total output of the two individuals is not additive, that is, 
together they produce in general more than separately. Each individual 
behaves in a way which he considers will maximize his satisfaction func- 
tion. Conditions are deduced for a certain relative equilibrium and for 
the stability of this equilibrium, i.e., conditions under which it will not 
“pay” the individual to decrease his efforts. In the absence of such condi- 
tions “exploitation” occurs which may or may not lead to total parasit- 
ism. Some forms of the inverse problem are considered, where the form 
of behavior is given and forms of the satisfaction function are deduced 
which lead to it. 


In a previous paper (Rapoport, 1947), hereafter referred to as 
I, we assumed that the total output of two individuals equals the sum 
of their individual outputs. As in I, we define x and y as the respec- 
tive efforts (and outputs) of two individuals X and Y, whose satis- 
faction functions are given by 


S;=R(@,y) + Hi(x,y), i=1,2, (1) 


where R;, essentially positive, represents the contribution to satis- 
faction by remuneration received, and H;, essentially negative, the 
detraction from satisfaction by effort expended. We shall call an opti- 
mal curve for X , the set of points (x, y) at which 0S,/ox = 0, and 
an optimal curve for Y, the set of points (~, y) at which 0S./oy = 0. 
Each individual then behaves in such a way (increases or decreases 
his output) as to reach a point on his optimal curve. 

We have shown in I that when total output is additive and is 
equally shared, the optimal curves for both individuals are represent- 
ed by the same straight line 


1 
oy -—2.-° (2) 
p 
Each point on that common optimal curve was shown to be un- 
stable in the sense that a small decrease of effort on the part of one in- 
41 
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dividual resulted in an increased effort on the part of the other ulti- 
mately leading to “complete parasitism” (see Figure 1). Furthermore, 


Y . 


a_— X+Y=|4-2 


Xx 


FIGURE 1 


Every point of the line x + y = 1/8 — 2 is an unstable “point of agreement.” 
If X decreases his output from E to EH’, Y will increase his output to #”. The 
point E” is a point of greater satisfaction for X , and he may repeat the process 
until his output is zero (complete parasitism). 


it was impossible under those conditions to find a point on the optimal 
curve which could be distinguished in any way as a “‘point of agree- 
ment” for the two individuals. : 

By a “point of agreement” we mean a point (%, ¥) which is a 
projection on the zy-plane of a point of the intersection curve of 
Si(a, y) and S,(a, y) which is characterized by certain maximum 
properties. If we demand that the “point of agreement” be a “total 
maximum” of both surfaces, to be defined below, the conditions on the 
partial derivatives of S; with respect to « and y are very stringent 
[see equations (54)]. If, on the other hand, we weaken these condi- 
tions and demand only that 


OS, OS» 07S, 07S» 
ou oy 0%? OY>” 


then every point of the optimal curve in Figure 1 is a “point of agree- 
ment.” 

We seek a situation in which the “point of agreement” satisfying 
equations (3) is unique at least in a suitable region. We shall then in- 
vestigate the conditions of stability at such points. 
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: In I we have considered six cases of motivation interactions of 
two individuals. We have numbered the cases in this paper consecu- 
tively with those of I. 

Case VII. The output is increased in joint production by a quan- 
tity proportional to the product of the outputs. Equal sharing. De- 
termination of output by simultaneous maximizing of S, and S.. 

The failure of the preceding system to give satisfactions greater 
than those achieved under individual production except in situations 
involving the exploitation of one individual by another, is traced, of 
course, to the decreased initiative on the part of the individuals to 
produce when maximum satisfactions are computed by partial differ- 
entiation (since each individual’s increased output is now split be- 
tween himself and his neighbor). In partial differentiation, no ac- 
count is taken of possible effects on the output of the second individual 
by the behavior of the first. We have seen how the introduction of 
contracts does account for such effects. We shall here consider incen- 
tives other than contracts. Let the efficiency of production increase 
with joint production (say by division of labor) by an amount pro- 
portional to the product of the individual efforts, so that for no effort 
on the part of one individual, the output reduces to the individual 
output of the other. Therefore, let the total output be given by 


T= Fa Yi OLY « (4) 
We shall refer to a as the cooperation coefficient. Then 
it ie i 
S,=log | 1 Yat Suid |-ee, 
(5) 
ety aige OU 
S2= log 1 ASE ly, 


The system 
Bi, Sg 


’ 


ox 7 oy 
: ES Pe ay. (6) 
an 
—_—____—_—_——— f=0; ———————_— $= 0 
at be Yt ary 2 oe + ¥ + xy 
now yields a unique pair of solutions for x and y, namely - 
a—26+Q 
Ly = : ? 7 
y 2aB (7) 


where 


Q= Vo?— 808? +48. 
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We choose the plus sign for the radical (the other solution is not a 
maximum) and substituting into equation (5) get 


a—2f+Q 
ae 5 


S,° = So =log(a + Q) — log 4 8? — (8) 


First let us note the limit of S;* as a approaches zero. That of the 


J Lee 
first two terms is log —-. Evaluating the last term by L’Hopital’s 
26 


Rule, we get 
a—2$8+Q 91 
gt Be (NB: lim Q=2 8). 


a0 


Lim 
a0 2 a 


1 1 
Hence Lim S;* = 2B ano + 6, which is exactly the value of S,* = S.”* 


in Case III of I for k = 4. But in Case III there was no way to come 
to this “agreement,” by simultaneous solution of the system 


0S; 0 Se 
On oy 


With the added motivation of the cooperation coefficient, no matter 
how small, such a possibility for agreement by simultaneous partial 
differentiation exists. 

Let us now see how S;* behaves with respect toa. We have 


a—4 pF a— 4 


eae CG) A o2 


ie B 
The limit of the first term as a > 0 is . The limit of the second 


term is indeterminate but can be evaluated by L’H6pital’s Rule. After 


; / awe 4 = sie 1 
two differentiations we find it to be — Hence 
him 28" 4 P88 +8 
m = 4 
eee vr >0 forp<#. 


Therefore, production is possible for 6 < 4 for all values of a > 0. 
It can also be shown that if Z is the maximizing value of the argu- 


0x 
ment, then ge > 0 for 6 < 4, that is, the cooperation coefficient is an 
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incentive for greater optimum output. 


: os 0S. 
The solutions of as = 0 and a = Oliealong two curves, name- 
Yy 


ly; 
1—2 6 — px; 1—28— pv. —a 22 
Pern Te a es 
afp2,+ B—ea ap t+ B 
An intersection of these two curves is the point (Z, y), Where 
a—2p+Q 
2a 


and Q = \/a?— 8a f? + 4 Bf. For the intersection to be real, we must 
have ‘ 


G2 6.0 f° AG? > 0. (9) 


That this is always the case for 6 < 4 is seen from the following con- 
siderations (NB: £ < 4 holds throughout the discussion, since it is a 
necessary condition for production to start). For 6 = i, the left side 
of relation (9) reduces to a2 —2 a+ 1, which is positive for all values 
of a. It remains to examine two cases: 1) 6 <4,a<4and2) p<}, 
a>4. Leta<4. Theno?—80/?+4#=o? —4 #(2a—1). But 
2a— 1 is negative if a < 4, hence Q? > 0. Now let a > $. Write 
Q? = ao? + 4 f2(1 — 2 a) and note that the second term is negative. 
But we have shown that Q? > 0 for 6 = 4. Reducing f reduces the 
only negative term, hence Q? > 0 if 6 < 4. In all cases considered 
the two optimal curves will intersect in a real point. Let us consider 
the “stability” of the system at that point. 

By stability of a system we understand a situation where it does 
not pay for an individual to lower his output. In this sense the sys- 
tem based on equal sharing without contract considered in Case III 
of I was not stable. For suppose X to be an “active” individual, i.e., 
capable of effecting changes in his behavior not merely on the basis 
of immediate increase of satisfaction, but with a view toward future 
events. Starting from any position of equilibrium in Case III, sup- 
pose X lowers his output. Since the position of equilibrium was a 
maximizing position for S,, this change has for its immediate result 
a decrease in S,, X’s satisfaction. But Y, acting entirely on the 
principle of immediate maximizing of S. , increases his output so as 
to get back to the optimal line (See Figure 1). Thus a new point of 
equilibrium E” is reached, at which point the satisfaction of X is 
greater than at FE, while that of Y is less. This process can theoreti- 
cally continue until X has become parasitic on Y, barring such phe- 
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nomena as the increase of “activeness” with decrease of satisfaction, 
which would cause Y ultimately to resist such a process. In a stable 
system, on the other hand, a decrease in output on the part of one 
individual would have for its end result a decrease in his satisfaction 
function. Nothing has been said here about the effects of increasing 
output on the part of an active individual. Theoretically, situations 
can exist where such an increase causes an increase in the output of 
the other, and so on indefinitely. Such cases are of no interest be- 
cause of the obvious physical limitations to increase of effort, limita- 
tions not expressed in the equations and therefore extraneous to the 
problem. 

In the present case, the situation is considerably different from 
that in Case III. We are dealing with two separate, nonlinear opti- 
mal curves, instead of two coincident linear ones (see Figure 2). 


axBX+B-a 


2 =o 2 4(a= 2x 
aBX+B 


FIGURE 2 


Optimal curves in the case of non-additive total i 

c n- output. If X decreases h 

Sane dea oe E to E’, Y will bring his own output to the point #” on hewn 
ptimal curve. For a > 3, this results in a decrease of satisfaction for X . 


if X now gets off the point of equilibrium by decreasing his ef- 
fort, Y will tend to bring the production back to his own optimal 
curve, distinct from that of X. To see how that will affect X’s satis- 
faction function ultimately, we must take the directional derivative 
of S, in the direction of the S, optimal curve. Let us compute the 
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slopes of the optimal curves at #, 7, their point of intersection. We 
have 


Pee ie 1) 
ein tia ee. 
62(2a—1) 
> 11 
see (Circ at 8) ~ 
4 B2?(2a—1) ; 
phere . : 
y Bs eo ahs (12) 
_ 46(2a—1) 2 
Y 2 aa (G0)? (13) 


Yo=y 
If 6, and 4, are the respective inclinations of the tangents to the 
optimal curves at (Z, 7), we have 
tan@=y',| _; 
W=y 
tan @.=y'2| _ 
YoY 
Also sin 6, 2 0 and sin 6, 2 0 since 6, and 62 are taken between 0 and 
a. 
We seek the directional derivative of S, in the direction 6: 


0S, 0S, 0S, ; 
— COS\U>.5- = SIO 
0s - One oy 
0S. 0S. lt+azx 
But a 0 at (“, ¥) and pt ———_—___—_————,, hence both de- 
Ox oy A ee a ty) es i 


rivatives are positive for all non-negative values of x and y. There- 


fore 
0S: 


0s 


Observe, however, that for a < 4, tan 6, <0. Therefore, fora <4, 
the 6. direction corresponds to decreasing x. First let a < 4. Then 
if x has decreased and the production has been restored to the opti- 
mal curve of Y, the point (x, y) has moved in the 6, direction. Since 
by equation (14) the directional derivative is always positive, S. has 
been increased by the decrease in +. Thus for a < }, the equilibrium 


IV 


0. (14) 
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is unstable at (@, 7). If, however, a > 4, tan 6 > 0, decreasing of x 
and restoration of (a, y) to the optimal curve of Y moves the point 
in the direction opposite to 6.. Since the directional derivative is still 
positive, this results in a decrease of satisfaction for X . He will there- 
fore tend not to decrease his effort. On the contrary, he may tend to 
increase it indefinitely, unless other factors (such as decline in effi- 
ciency for very high values of « and y, or increase of reluctance or 
any other stabilization of the system) interfere. We therefore state 

Theorem 7. If the satisfactions are given by equations (5), an 
optimum output can be agreed upon by simultaneous maxinuzing of 
the satisfaction functions partially with respect to the respective vari- 
ables. The system is stable at the optimum output if and only ifa 24. 

Case VIII. Determining the output from the point of view of X. 
Now let X be active, i.e., X considers the situation in advance and 
has some knowledge about the behavior of Y. Then X considers that 
Y will try to maximize his satisfaction function, i.e., try to bring it 
onto his own optimal curve. Therefore X takes the value of y given by 
0 Se : 
= — = (), thatas, 
oy 

Le bo oe 


eae api+tp 


Substituting this value into S, and differentiating now totally with re- 
spect to x gives for optimum output 


(15) 


we (16) 
The individual X lets Y determine his optimum output according to 
equation (15). Thus 

0? == a Bb 2.0 Bie eB 
= a : (17) 
Substituting the values of x and y given by expressions (16) and (17) 
into S,(“, y), we get 


21 


S,=log ( = pa 
S.=log (== JP E (18) 


Elementary computations show that, in effect, a stable equilibrium with 
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respect to X’s behavior has been reached. This point of stable equilib- 
rium could have been computed also from the preceding case by set- 
ting the proper directional derivative equal to zero. We can now ex- 
amine in detail the state of affairs determined by the “active indi- 
vidual” simply through his superior knowledge of the situation. In 
particular, we can compare the resulting satisfactions with those 
characteristic of the simultaneous solution equilibrium. We shall see 
what changes have occurred in Y’s (the passive individual’s) satis- 
faction function. Finally, we can compute the value of a for which 
it becomes profitable for one or both individuals to enter joint pro- 
duction under sharing. 

We note first of all from equation (16) thatifa=$,«=0.Since 
only non-negative values are considered, we see that the condition 
leads to parasitism by active X upon passive Y if a S 8. In that case, 


S, = log ( =) : 
S.=log ( = )-1 420, 


exactly as in the corresponding degenerate case (complete parasit- 
ism)in Case III. 

If a > 6, a point of equilibrium short of complete parasitism is 
reached. 

Theorem 8. Parasitism will occur in Case VIII if and only if 
a = p. 


Therefore, let a > 6. We have 


1 
0a ist hee Yor a 
og 1 
da a = 


and we see that for greater a, X will choose a point of greater effort 
for himself, but his resulting satisfaction will also be greater. 

Let us now compare S, as given by expression (18) with that 
given by equation (8) and see under what conditions the former will 
exceed the latter; in other words, we wish to see whether the fore- 
sight exhibited by X will always pay under these conditions. 

The difference in the two satisfactions is 
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a a—~p a—26+Q 
— es Se ogg h? ee 
log op . log(a + Q) + log 46 oa 
1 
= log 8-24 2 =D, (a) 
aie 


Note that D,(4) = 0 since Q(4) =4=a. 
Consider now that 


eS) 2 SO) 
ONO 8 Qa + a) 


Examining Q(a) = \/a? — 4 f?(2 a — 1), we see that a < Q if, and 
only if, a <4. Hence 


O@) 


FIGURE 3 


Comparison of advantages to X and Y, when X ef ip” 
conditions of Case VIL g s n X assumes “leadership” under 
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D,'(a) <0 for a<4, 

D,' (a) =0 fora=}, 

D(a) >0 for $<a<1, 

D(a) <0 for a>1, 
Lim D'(a) =0. 

This examination of D,'(a), plus the fact that D(4) = 0, enables 
us at once to draw a graph of D,(a) (see Figure 3). We see that the 
satisfaction of X arrived at by his computation (Case VIII) is great- 
er than that arrived at by simultaneous partial differentiation (Case 
VII). However, referring to Theorem 7, we see that for a < 1/2, 
A tends to decrease his efforts to reach the new equilibrium, while 
for a > 1/2, he tends to increase them. In the first case, he will ar- 
rive at greater satisfaction by simply exploiting Y; in the second 
case, he himself will also increase his efforts. 

It remains to consider how Y’s satisfaction will be affected by 
the manipulations of X. Again as above we take the difference in 
Y’s satisfactions as given by equations (8) and (18). The difference 
is 


Q 1 2 B? B? 
Dz (a) =log(2 a) —log(a + Q) +——-—=+ eae 
LS 4 a a? 
(19) 
§ (2.01) 
=D;(c) +——_—. 
a 


The superimposed graphs of D,(a) and D,(a) are shown in Fig- 
ure 38. Under the conditions of Case VIII, it is always advantageous 
for X to assume “leadership”. The degree of advantage is given by 
the value of D,(a). The value of D.(a) represents the advantage to 
Y in accepting the leadership of X . This becomes positive for a > 4. 
When « > 4, it is even more advantageous for Y to accept X’s leader- 
ship than it is for X to assume it. 

It remains to consider under what conditions X will join the sys- 
tem in Case VII and in Case VIII. When X works alone, his maximal 
satisfaction is given by equation (5) of Case I of I. 


Sax =—logp+f—1. 
This subtracted from X’s satisfaction gives in Case VIII [cf. equa- 
tion (8) ] 
a—28+Q 


D,(a, 8) =log(a + Q) aN a cea gees a ail big ek at (20) 
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Thus D,(a, 8) > Oif 
at+Q _ 2ap+Q—2p—a (21) 
oe ( 4B pp 20 


For the critical value a = 1/2, the condition on f is especially simple, 
namely, 


log 4 p<. (22) 


Since the satisfactions in Case VII are equal, similar conditions 
hold for Y. We summarize these results in 

Theorem 10. Sufficient motivation for either individual to co- 
operate and share under the conditions of Case VII is given by rela- 
tion (21). For the minimum value of a which gives a stable equilrb- 
rium (a= 1/2), the upper limit on B is given by relation (22). 

For Case VIII, the conditions are less stringent for X , since, as 
we have seen, X always increases his satisfaction over that of Case 
VII. We have for X 


p 
DD 9 ] ee (28) 
(a, 6)'= ary. F B. 
This is positive if 
Qa Op as 
joz ( = )> P uh (24) 
26 a 
oD Leet : 
aaNet tage nay ae 7 <0 fora > 8. Hence if-we let £ 
Qa 


be 1/2, its maximum value, we can deduce a simple sufficient condi- 
tion on a for the non-parasitic case, namely 


1 
log a + — > 1/2. (25) 
2a 


We see that if 8 is sufficiently small, relation (24) will always hold 
if a is fixed and finite, but not necessarily if a also approaches zero. 


aie a 
In the parasitic case S, = log ( De ) : 


D,(a, ) =Iog ( § 2 )+1—10g2— B. (26) 


Ifa= 8, the arrangement will pay if B < 1— log 2, which i is 
the condition on f in the parasitic situation in Case III of L 

We should expect the condition for Y to join X to be more strin- 
gent in Case VIII since Y is to play a passive role. However, this is 
not always the case. We have for Y , 
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Oo ps2 a foe B* 
2 


a 


D. (a, 6) = log ( = ) + log B— 6 +1. (27) 


This is positive if 
a = 2 
jog ( =~ \> SFP 1—2a), 28 
a Fai of hae | (28) 
and we have the curious result that for a > 1 /2 it is even more ad- 
vantageous for Y to join X under X’s leadership than it is for X to 
enter the arrangement (see Figure 3). 

Case IX. The individuals increase their output by imitation. Let 
us now consider a situation similar to that in Case III of I, except 
that in addition to the effort which each individual expends to satisfy 
his own satisfaction function, he also expends additional effort in pro- 


portion to the effort of the other. In other words, we are dealing with 
imitative behavior. The equations connecting the variables are now 


1 
L=-—2—y+ky, 

B 

t (29) 
f= — 2 at he. 


The optimal curves are now two distinct straight lines intersecting at 
1—2, 
£= y¥ = ———— > 0 (30) 
B(2—k) 
if k <2. The situation has no meaning for k > 2. 
If the point of intersection of the optimal curves determines the 
efforts of X and Y, we have 


ths yong ) 1—2 8 
HS, = 1 (=. rar , (31) 
See = g(9 =) ) = i 
which for k = 0 reduces to S* of Case III of I under equal division of 


labor. Thus the additional motivation of imitation is, like the coopera- 
tion coefficient, a “stabilizing” influence in that a common optimum 


can be uniquely found. 
Comparing the S* of this case with Snax of Case I, we see that 
the imitation coefficient is sufficient motivation for cooperation if 
ue ) 1—28 
B(2—k) 2—k 


D(B, k) = S* — Sou = log ( 


+ log 6—p +1 0, 
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that is, if 
jor (3) 4 IE oo. (32) 
2—k 2—k 


Let us examine the situation for stability and for conditions for 
complete parasitism. If X takes the initiative, he computes the be- 
havior of Y from the second of the equations (29). Then 


ee eae \ es 
$,=108 (— Bu, (33) 
k—1 
_—=—_. 34) 
x ER ( 


The condition for stability of the point (@, y) determined in equa- 
tions (30) is, therefore, 


k—-1, 1-28 - 


k Bx B(I2—k)” 
(2—k)-(k—1) 2k—2 pk; 
(= WAS2 6h). (35) 
The region of stability is shown in Figure 4. From equation (34) 
B 


pant = V2 Te 


The shaded area represents the region of stability for Case IX. For the 
points (k, @) in the region assumption of “leadership” by X will not result in 
exploitation. The range of @ is largest at k = V2. For smaller values of k Xe 
will tend to exploit Y because he cannot make Y work sufficiently merely by “set- 
ting an example”. For larger values of k, his own tendency to imitate Y limits 
the tendency to exploit. 

-B. In spite of imitative behavior in Case IX, the satisfaction function is 
assumed to be of the same form as in the preceding cases. 
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we see that parasitism occurs for k = 1. From inequality (35) we 
we see that stability is insured if 6 is sufficiently large (provided 
k > 1, since for k < 1, inequality (35) can never hold, B being < 4). 
Hence 

; Theorem 11. Sufficient motive for cooperation under the econdi- 
tions of Case IX is given by expressions (32). Parasitism does not 
occur fork > 1, and the equilibrium is stable if fh ws suf ficiently large 
to satisfy inequality (35). 

Remark on Case IX. This case is peculiar in that conditions for 
stability are exhibited as conditions on 6 (provided parasitism does 
not occur), where a large value of f insures stability. An interpreta- 
tion of this phenomenon can be given roughly as follows: If £ is 
large, Y will not increase his production sufficiently when X puts the 
pressure on him by lowering his own production. It will pay X rath- 
er to depend on Y’s imitation coefficient and have him produce by 
“setting an example” rather than by pressure. 


Now we will discuss an inverse problem. 


Case X. Forms of satisfaction functions, from which equations 
(29) can be deduced. 

In the preceding case, instead of specifying a satisfaction func- 
tion, we have specified the final behavior. This may be considered as 
a situation defining an inverse problem. For what forms of the satis- 
faction function will the given behavior be a maximizing behavior? 
A completely general solution is probably not interesting. It would 
be analogous to the problem of finding all functions of x which have 
a maximum at « = £. But with proper restrictions the problem does 
yield interesting results. 

We shall state the problem in the following form. Let S, and S, 
be given by equation (1), where FR; is a logarithmic function of a 
function of x and y (R; = log [f(x, y)]) and E; is a polynomial. 
What change in S; from its form given by equations (5) will result 
in a maximizing behavior described by equation (29)? For simplicity 
we leta=0O. 

Let us suppose this change is produced by two functions 2, (x, y) 
and z.(x , y) somehow interacting with the S; of equations (5). Even 
if the equations are to preserve their logarithmic-polynomial form, 
-2(x, y) can interact with the S; of equations (5) in many different 
ways. It can enter as an additive term, a factor, or an exponential 
in the argument of the logarithm, etc., or as an additive term or a fac- 
tor in E , etc. Depending on the nature of this interaction, we are led to 
innumerable forms of the satisfaction function, all of them leading to 
equations (29). While it is true, therefore, that a particular form of 
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the satisfaction function determines a particular form of maximiz- 
ing behavior, the converse is by no means true. If at any time prob- 
lems dealing with actual behavior are formulated on the basis of a 
maximizing behavior, it is, of course, the behavior which will be ob- 
served, and so the mathematical sociologist will be confronted with 
the ambiguous inverse problem instead of the straightforward direct 
one. As we shall show, however, the case is by no means hopeless. 
Other criteria besides the maximizing behavior enable us to accept 
certain of the possible forms of the satisfaction function and to re- 
ject others. 
We first suppose that the change is an additive one in E: 


Lt ¥ 
S,=log ( 1 ha aes )-s2+ale,y) : 
cia (26) 
S.=log (1 ore —By+“(x,y). 
Then 
rs) 4 itt 0 41 
= ———_ — f+ a= ()e 
se ae ail wesc 0x 
gee pene 2+kh ah 
= -— — = = -—-— y— z 
; pax y Y 
The general solution of the partial differential equation (37 \ris 
2, = pkay + fly), (38) 
where f(y) is an arbitrary function. 
Similarly, 
2 = Bkay + g(x), (39) 


with g(x) an arbitrary function. In the simplest case fw=2(%)= 
0, 2% = 2 = pkxy. The form of S; is somewhat like that in Cases 
VII and VIII except that the xy term is added to E instead of to the 
argument of R. We have 


te 
S,=log (1 +=) = bet ples 


(40) 
x 
S2 = log (1 a2 a ) — by + pkay. 
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Following the methods of the preceding sections, we could de- 
duce analogous conditions of stability, sufficient motivation, etc. We 
shall not do this here, however, but instead inquire into the form of 
& and z, if they are to enter as additive terms in the argument of 
Rk. (Note that if we consider z as a multiplicative factor in the ar- 
gument of # , we obtain essentially the same results as in the preced- 


ing case.) We shall for convenience consider z/2 as the additive term. 
Then 


fae) a 
S,=log (1 + ———__— )= pe: 
0 z 
as hi a 
x 
“= Seen memes dae 
OL 2 oy ey 
(41) 
E 9 ioe ae ae 
=-—2—y-—z _ Sac OP : 
ety dans i fiP %.5 6 ) 
— — pa, = phy. (42) 
Cx 
The general solution of the partial differential equation (42) is 
2 = e8 f,(y) — ky, (43) 
where f,(y) is an arbitrary function. Similarly, 
Z — ef fo{a) — ke . (44) 


In spite of the restrictions we have imposed on 2, and 22, we see that 
the solutions, being solutions of partial differential equations, are 
characterized by great generality. But we can further restrict the 
character of these solutions by additional considerations. For exam- 
ple, it is reasonable to demand that z, increase monotonically with y. 
Hence 


2 = eb f',(y) —k>0 (45) 
OY 


for all values of x. Then, since only non-negative values of x are 


considered, 
f.a)e2k.. 


Or we may weaken this restriction and demand that merely the part 
of the argument of R containing y should increase monotonically with 


y. Then 
f(y) ie bs 
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In the parasitic case f'1(y) can actually be negative (k <1). Taking 
a simple function for f,, f:(y) = ky, we have 


2,—= ky (e6 —1) 20 (46) 


for all values of x and y. We might further inquire whether there 
exist forms of the functions f,(y) and f2(y) which would make 
Z, = 2, i.e., 2; symmetric in x and y. If such a form could be found, 
we should be dealing with a situation where the symmetric character 
of R; is preserved even with the introduction of imitative behavior. 
It appears, however, that the functional equation z,(x , y) = #(",y), 
where z, and z, are of the form given by equations (43) and (44), 
has no solutions except for k = 0. Hence 

Theorem 11. If imitative behavior given by equations (29) is 
brought about by an additive term in the argument of R, then R ts 
not symmetric in x andy. : 

Case XI. Each individual desires a maximum output of the other 
over a period of time. 

Introducing arbitrary functions into the satisfaction function, as 
we have done in the inverse problem, makes it possible to impose a 
large variety of additional-conditions. The forms of the satisfaction 
function we have considered in the direct problem were such that the 
equations 


0 Sy 0 Se 
—=0 and =0 
oY. x 


had no solutions. The psychological interpretation of this is, of course, 
that it was always in the interest of one individual to have the other 
one work as much as possible. The negative “effort” term of X de- 
pended only on «. We could, of course, introduce many other forms 
of the satisfaction function which would have the effort term depend 


aS 
also on y and thus make possible the solution of ae = 0. But which 
y 


of the infinity of such satisfaction functions shall we choose? 


Abstract and formal as our treatment is, we can try to keep at 
least a semblance of reality by deducing the forms of the satisfaction 
function from situations which describe an individual’s interest or 
behavior. This we did in the inverse problem. We started with be- 
havior (e.g., observed effect of increasing productivity of one indi- 
vidual on that of another) and deduced from it satisfaction functions 
determined by various forms of z. If we can consider the reluctance 
not a constant as we have hitherto considered it, but as a variable 
which is monotonically decreasing with respect to y for X and with 
respect to x for Y , then we have, for example, 
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bs = Bo(1— ky), 

Bs = Bo (1 — ka). 
Equations (40) then become a feasible interpretation of the observed 
behavior described by equation (29). 

Let us now endow our individuals with a foresight. Their desires 
are now not so much for the maximum immediate rate of production 
of the other individual but for a maximum output over a period of 
time. 

To see how such “foresight” operates, consider that X now pos- 
sesses a work animal or a machine Y , which produces at a variable 
rate y = y(t). This rate is monotonically decreasing with respect to 
time (because the machine is being worn out) and the rate of its de- 
crease depends on the initial effort y, and on the variable effort y . This 


can be described by a differential equation with a boundary condition, 
for example, 


(47) 


y=—yy; y(0)=%. (48) 
The rate of output y is now a function of ¢. The solution of the equa- 
tion (48) is 


Y=Yo ev", (49) 
In a unit of time 0 = t = 1, Y will yield the total output 
. 1 
w= if ydt =— (1—e“») = W (yw). (50) 
\ 0 Yo 
Yo 


FIGURE 5: 
Graphical solution of the equation e-¥’(i + 2y,?) =1. 
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If X wishes to maximize W with respect to the initial value of y, 
then differentiating with respect'to y, we get 


dW 
dyo 
e-vo(1 + 2y%) =1. (52) 


An obvious root of equation (52) is yo = 0. We see, however, that 
this is an inflection point at W = 0. Another root of equation (52) 
exists as is readily seen by plotting 


1 — evo — 2y?, ew = 0; (51) 


Uw—e% and w%———. 
3 1 + 2y%, 


(See Figure 5.) Note that for yo = 0, w= %=1. For small posi- 
tive values of Yo, Us > Us Since 


du, ‘due 


= =0 at y—0, 
dy. Yo : 
and 
Cu, a 
: bs at y= 0. 
dy?) — dy? 


On the other hand, for large values of yo, wu: < uw. The root of 
Yo of equation (52) can be calculated from the value of wu, where 
WU = uU,. This reot can be shown to give a maximum W. We shall 
therefore call it y¥. Similarly, we can compute #, which by the sym- 
metry of the situation equals 7. 


We can now state the inverse problem as follows: What forms of 
the satisfaction function will insure 


0 Sa 
oy 
0 So 
Ox 


(53) 


w=t 


Finally, we.can ask for the forms of the satisfaction function 
which will make the wishes of each individual with respect to his own 
and his neighbor’s optimum output coincide. This is tantamount to 
the following conditions at the point (#, 9): 
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0S, 0S, 08, 08, =) 


Oy oy Ox oy 
eS. \*_ (aS, 2S, 
-|}<0, 
oxoy 02? Oy? 
Ree \s se / 08 8, 08.8: 
Oy 
ox oy rae wD 
he Si 0? Ds 0 So 0 S» 


0, ar Os <0. 
0 x? 0 y? 0 20? dy? 


’ 


(54) 


<0: 


A society of two individuals characterized by satisfaction func- 
tions which satisfy the conditions (54) would be absolutely stable in 
the sense that it would not be of interest to either individual to change 
either the rate of his own output or the rate of the other’s output. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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Using expressions derived in previous papers, the author investigates 
the behavior of a cell immersed in an infinite medium, under the influ- 
ence of diffusion of a single solute and flow of water. The effect of hy- 
drostatic pressure on the system is taken into account. It is found that, 
depending on the values of certain parameters, the cell can collapse, 
burst, reach a stationary stable state, or execute undamped oscillations; 
a cell must burst or collapse unless its volume is an increasing function 
of internal pressure, and it can execute stable oscillations only if its 
membrane acts as a “potential well” to the molecules of the solute. 


In previous publications (Bloch, 1944, 1946), equations were de- 
rived giving the permeability h (to flow of a solute) of a membrane 
separating two regions 0 and 1. Equation (22) of the latter paper is 


J D 2B 3/2 —z u ev-% — y er? 
h=—— == (1— >) e= ) ca 
€o—C, 2b Bilis u—v er” — ev-u 


The symbols are all defined in the papers cited, except for ¢., which 
here replaces the V, used previously. In particular, J is defined as 
the flow of the solute, in grams per square centimeter of membrane 
per second, from region 0 into region 1; the concentrations of the 
solute, c. and c, , and the hydrostatic pressures, P, and P,, in regions 
0 and 1 respectively, enter the equation through the quantities u,v, 
y , and z, defined below: 


2bH _ 2bK 
“u= Co» ae Ps 
Dy 0 
| (2) 
2bH _ 20K 
p= Gy, A ‘cake 
0 0 


Here b and D, are positive constants; H and K, also positive con- 
stants, are defined in terms of J, the flow of water in grams per 
square centimeter of the membrane per second, from region 0 into 
region 1: 
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I=H(¢,— Co) + K(Po— P,) 


(3) 
=> [(v—u) + (y—2)]. 


Equation (1) includes the effect on the motion of the solute of water 
flow due to both hydrostatic and osmotic pressure gradients; how- 
ever, this equation has so far been applied only in the special case 
of equal hydrostatic pressures on the two sides of the membrane 
(y =z). The purpose of the present paper is to investigate the more 
general case in which y # 2. 

The simplest system which seems likely to be interesting is a 
single cell in an infinite medium. If the interior of the cell is taken 
to be region 1, the equations of state of the system are 

u = const. 
and (4) 
- y = const. 


In addition, the following equations will be needed: 


2bH 
(Vv)'=AI =Ah(u—v), » 4) 
Do 
and 
DA 
VS Al x (v—uty—z). (6) 


Here V is the volume in the cell occupied by the solution, A is the 
surface area of the cell, and a prime denotes differentiation with re- 
spect to time. Equation (5) states that the rate of change of the 
amount of solute in the cell is equal to the rate at which solute enters 
the cell through the membrane; equation (6) states that the rate of 
increase of the volume of the cell is equal to the rate at which water 
enters it. Obviously equation (5) is not accurate if there is metabo- 
lism of the solute in the cell; however, the other metabolites are likely 
to have somewhat the same effect on water flow as the solute which 
is considered, and explicit treatment of more than one solute compli- 
cates the problem greatly. Equation (6) should be reasonably accu- 
rate unless the concentrations of solutes are extremely high. It will 
be assumed that V depends only on z — y, or in this case only on z; 
ae z' has the same sign as V’ if dV/dz > 0, opposite sign if dV/dz < 

It will now be possible to find the equations of the curves v' = 0 
and 2’ = 0 in the vz-plane defined by equations (4). If v’ = 0 , equa- 
tions (5) and (6) combine to give 
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Do 
pe YOU hee), (7) 


or, from equation (1), v’ = 0 if 


U 
1— fp- 
A a em de (8) 
where 
p=(1 aie J" >0 9 
SkT, : (9) 


\ 

Two cases can be distinguished: 6 < 1 and 6B > 1. These cases corre- 
spond, respectively, to the membrane’s being a potential barrier and 
a potential well. When § < 1, equation (8) shows that v’ cannot 
vanish for real values of z unless v 2 fu; v' vanishes along a curve 
of positive slope which is asymptotic to the line v = fu as z approaches 
— co and whose slope dz/dv approaches + 1 as v and z approach + o. 
When £ > 1, the curve v’ = 0 lies between the lines v = 0 and v = fu; 
as v approaches zero, z approaches + o , while as v approaches pu, 
z approaches — oo. 

The curve z’ = 0 can be obtained on the assumption that V de- 
pends only on z and that dz/dV does not vanish for any realizable 
value of z. Thus z’= 0 when V’ = 0, or, by equation (6), 2 = 0 if 


zg=v—uty. (10) 


This equation defines a straight line of unit slope in the vz-plane. 


Z’=0 
z=z* 


' 
' 
' 
‘ 
iS z=z7** 


1 

1] 

' 

! 

uw \ Vv 
FIGURE 1 
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B >i, dV/dz>0 


FIGURE 2 


Figures 1 and 2 show graphs of the curves v’ = 0 and z’ = 0 in 
the first quadrant of the vz-plane (negative values of v or z are not 
possible). Figure 1 depicts the situation for 6 < 1; Figure 2, that for 
f>1. For definiteness wu is taken greater than y in both cases; there 
seems no reason to believe that the relation between uw and y has any 
qualitative effect on the behavior of the system. In all cases the two 
curves v’ = 0 and 2’ = 0 intersect at the point v = u, z= y and no- 
where else. Clearly when there is no flow of water (z’ = 0), v’ is 
positive when v < x and is negative when v > u; thus v’ > 0 to the 
left and v’ < 0 to the right of the curve v’ = 0. Also, when v = u, 
V’ is negative or positive according as z is greater or less than y; 
z > 0 to the right and z’ < 0 to the left of the line 2’ = 0 if dV/dz > 0, 
has opposite signs if dV/dz < 0. 

The first case to be considered is that in which V is an increasing 
function of z. In that case, if a representative point is on the line 
z’ = 0 in either figure, its motion must be parallel to the v-axis, to- 
ward the right if the point is below and toward the left if the point 
is above the intersection of the two curves. If a point is on the curve 
v’' = 0, its motion is parallel_to the z-axis, upward if the point is 
below and downward if the point is above the line z’ = 0. If a repre- 
sentative point is at the intersection of the curves v’ = 0 and 2’ = 0 
(Le, atv =u, z= y), it remains there; this is the only possible sta- 
tionary state for the system. A point at the right of both curves moves 
upward and toward the left; a point at the left of both curves moves 
downward and toward the right; a point between the two curves 
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moves upward and toward the right if it is below the line 2’ = 0 - 
downward and toward the left if it is above that line. The situations 
ees above are represented by the small arrows in Figures 1 
and 2. 

In the case being considered (dV/dz > 0), when Bb <1, a repre- 
sentative point in Figure 1 which is at any time between the curves 
v’ = 0 and z’ = 0 remains between these curves and moves toward 
their intersection, which represents a stable configuration. A repre- 
sentative point which starts at the right or the left of the two curves 
moves toward the curves along a path of negative slope. Unless a 
catastrophe occurs first, the point eventually crosses one of the curves 
and goes to the stable point (wu, y). However, there are two limiting 
values of hydrostatic pressure within the cell, hence two values of z: 
one value, z = z*, corresponds to rupture of the cell membrane, while 
the other value, z = z**, corresponds to collapse of the cell; or, if the 
cell does not collapse, z** may correspond to the vapor pressure of the 
solution in the cell. Thus, in the case depicted, if the initial state of 
the system is represented by a point at the right of the two curves 
in Figure 1, this point ultimately reaches stability at (u, y) or else 
crosses the line z = 2* so that the cell membrane bursts. If the repre- 
sentative point starts at the left of the two curves, it reaches the 
stable position (uw, y), or else crosses the line z = z**, so that a non- 
rigid cell collapses; if the cell is rigid, it may be that the internal 
pressure falls to the vapor-pressure of the solution in the cell (unless 
the constants of the membrane are changed first) , so that some water- 
vapor is formed in the top of the cell, and water flows in the top and 
out the bottom keeping the internal concentration of solute at some 
constant value. 

When £ > 1, the situation is more complicated. Figure 2 shows 
that a representative point tends to circulate around the equilibrium 
‘position (u, y); however, it is not clear from the Figure whether 
such a point moves toward or away from (uw, y) or, if it does approach 
equilibrium, how many oscillations the system executes between any 
given initial state and the final state. Complete answers to these ques- 
tions could be found only from the solution of a differential equation 
which will be derived later. 

However, it is possible to show by a simple argument that in 
some cases the system does not “blow up”. Imagine a vertical straight 
line-segment in Figure 2, drawn upward from a point (v1 > 21) at the 
right of (w, y) on the curve v’ = 0, until it intersects the line 2’ = 0; 
any representative point which is initially in the region bounded by 
this segment and the curves v’ = 0 and z = 0 can leave this region 
only by crossing the line 2’ = 0. Now suppose there is a line segment 
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drawn from the intersection of the first segment and the line 2’ = 0 
horizontally toward the left until it intersects the curve v = 0 at 
some point (v2, 22) at the left of (wu, y); any representative point 
which is initially in the region bounded by this second line segment 
and the curves v’ = 0 and z’ = 0 can leave this region only by cross- 
ing the curve v' = 0. The construction of line segments can be con- 
tinued: a third segment can be drawn downward from (?, Zo) 
until it intersects the line 2’ = 0 at the point (v2, 23), and a fourth 
segment can be constructed toward the right from (v2, Zs) until it 
intersects the curve v’' = 0 at the point (v;, 2s). It is clear that no 
representative point can cross any of these lines in the outward di- 
rection; if the point (v,, 2,) is not between (wu, y) and (v3, 23) on 
the curve v' = 0, a representative point that starts in the region 
bounded by the four line segments and the portion of the curve v' = 0 
that lies between (v,, 2:) and (v;, 2;) can never leave that region. 
If none of the four segments intersects either of the lines z = z* and 
z = 2**, any system whose representative point is at any time in the 
region described above must either approach stability at (w, y) or 
eventually execute undamped oscillations about that point. If (v1, 21) 
is between (wu, y) and (v;, 23) it is not possible to conclude from this 
type of reasoning that the representative point is confined to the re- 
gion in question. : ' 
From the geometry of the figure it is evident that (v;, z;) coin- 
cides with (v,, 2:) when (2, — 2:)/(v2 — v1) =—1, and is between 
(w,y) and (v1, 21) when (2. — 2) /(V2— v1) <—1. The curve v' = 0 
lies entirely between the z-axis and the line v = fu; its slope approach- 
es — oo as z approaches either + o or — o; therefore, if (v,, 2) is 
far enough down the curve (perhaps below the v-axis), the average 
slope of the curve between (v,, 2) and any other point on the curve 
is less than — 1. Thus there is in every case a finite region of the 
ve-plane which a representative point can never leave; obviously, 
however, if the smallest such region contains parts of the lines z = z* 
and z = z™, the cell in question may still burst or collapse. In some 
cases, viz., when 
3 
2 


the slope of the curve v' = 0 is everywhere less than —1, so that the 
construction of line segments described above will, if continued, pro- 
duce a spiral curve which has (u, y) as a limit point. No represen- 
tative point can cross this spiral in the outward direction; therefore, 
in this case any representative point that starts in such a position 
that it does not cross either of the lines zg = z* and z = z** will ap- 


VBu< (11) 
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proach stability at (wu, y), though it may revolve about (u,y) sev- 
eral times in the process. Clearly the condition (11) is too strong; 
even if condition (11) does not hold, it may be that any representa- 
tive point that does not cross z = z* or z = z** will approach (u, y). 
The other alternative is that its path approach some stable closed 
curve from either inside or outside that curve, so that the motion of 
the cell approaches oscillation of constant amplitude. ; 
Some further idea of the behavior of the system may be gained 
from the following analysis. Equations (1), (5), and (6) imply that 


DA (z uU ev" — y er” 
v= 2 wv —y +2) v2 I], (12) 
2bV (z) ev-4 — @2-v 
and 
DA (2) 
a Se Vey anne ) ns 13 
2! = a (UU y+ 2) (13) 
2b — 
dz 
whence 
a= —| |. (14) 
dz Vdz a 


Equation (14) determines the path of any representative point in 
‘terms of its initial position. However, the exact solution would be 
quite hard to obtain, and would very likely be too complicated to be 
of use. Therefore, the procedure adopted will be to reduce equations 
(12) and (13) to linear equations, which are approximately valid in 
the neighborhood of (vu, y). The functions v’ and z’ can be expanded 
in Taylor series about the point (uw, y): 


yu (v,z) =v (u,y) + @—w(— ) + c—»( J tes 


ad ms te (15) 
Oe 2) 2 (sy). + (v—w( — ) us @—w( >> eee 


These two series both converge to the functions they are being used 
to represent over a region of which (wu, y) is an interior point. There- 
fore, if each series is cut off after the terms which are linear in 
(v —u) and (z—y), the values of (v — u) and (z— y) can be taken 
small enough so the series approximate the appropriate functions to 
any required degree of accuracy. Thus, since v'(u, yy =2(u,y)=0, 
the equations 


Uy 
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dv av’ ou’ 
aia (0 ae ok @ ae 
Legere (v—u) ( Ae \e ( w(s a ke 
and (16) 


dz 02 02 
oe —-— = (y—Uu ——. +@—w( ) 
° dt ( )( Ov Je 02 UY 


can be used to investigate whether (wu, y) is a position of stable or 
unstable equilibrium. When the derivatives in these equations are 
evaluated, the equations become 


dv DA (y) le 


and (17) 


The solutions of these equations are 


VY = Ay, et '+ dy. er" 
and ; (18) 
2 Oy Cae 
where the a;; are constants whose values are irrelevant to the present 
discussion, J, and 4, are the two quantities 


zee 2 Ip [w(up—u—p)—1+ 


Viv p—u—B) +1 —Fy@up—w | 


and Y is the logarithmic derivative of V(z) when z= y. 

“If the radicand in equation (19) is negative, the v and z coor- 
dinates of a representative point oscillate about and y respectively ; 
the amplitude of the oscillation increases or decreases with increas- 
ing time according as the part of 4, or 4 outside the radical is posi- 
tive or negative. When f > 1, the radicand is negative if 


Ai, 42 = 
(19) 


1 eS T as 
AMER ones rene TS SS (Vp Pr 1)*5* 720) 


if 6 < 1, the radicand is positive for all positive wu, or oscillation of 
the cell can not occur, as was indicated by Figure 1. If dV/dz <0, 
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y < 0, and there can be no oscillation; this case will be discussed 
later. 

The quantity in front of the radical in equation (19) is positive 
if y > 0 and 


pepe 1 WF vt vit 
vB-1) yB—1 2 


it is negative if y > 0 and the inequality sign is reversed. Thus if 
expressions (20) and (21) hold, the motion of the cell approaches un- 
damped oscillation unless the representative point crosses one of the 
lines z = 2* and z = 2**; if expression (20) but not (21) holds, the cell 
oscillates toward equilibrium with v = u and z=y. The period of the 
oscillations in the neighborhood of (uw, y) is approximately restricted 
by 

/ 


(21) 


4nbV(y) Vy 4nbyr y 
DA(y)VB Do pe 


where r is one of the linear dimensions of the cell when internal and 
external pressures are equal, and » is the dimensionless ratio of V (y) 
to A(y)r, assumed independent of y. In the case of a spherical cell, 
if r is the radius of the cell when internal and external pressures are 
equal, y= +4. The period of oscillation 7 has its minimum value, giv- 
en by expression (22) with the equality sign, when 

wht 


i; 
y(B—1) 


7 increases as wu increases or decreases from this value, approaching 
infinity as u approaches either of the limits given by condition (20). 
Of course, expression (22) or any other relation derived from equa- 
tions (16) can give at best only a very rough idea of the period of 
any oscillation that may occur far from the point (wu, y). 
When £ < 1 or wy < 0, or when condition (20) is not satisfied, 
the solutions of equations (16) are not oscillatory, or the motion of 
a representative point is not oscillatory in the immediate neighbor- 
hood of (uw, y) where equations (16) are approximately valid. How- 
ever, if (wu, y) is not a position of stable equilibrium, that is, when 
inequality (21) is satisfied, and 8 > 1 and y > 0, the straight line 
construction described above shows that the representative point has 
a stable closed path around (uw, y) on which its motion is escillatory, 
whether or not there is oscillation in the immediate neighborhood of 
(u, y). Thus expression (21) is the condition for the existence of 


IIV 


- 


(22) 
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stable oscillations of the cell, and when condition (21) is not satis- 
fied, 6 being greater than unity and dV/dz being positive, the cell 
approaches a stationary state in which hydrostatic pressure and so- 
lute concentration are each the same inside as outside. Incidentally, 
as is to be expected, the very strong condition (11) implies that weak- 
er condition (21) does not hold. It will be recalled that when 6 < 1 
and dV/dz > 0 the cell always approaches the stationary state. The 
foregoing remarks apply, of course, only if in the process of approach- 
ing the final state the cell does not burst or collapse. 


B >',.dV/dz<o 


FIGURE 4 


INGRAM BLOCH 13 


When dV /dz < 0, the situation is as depicted in Figure 3 (8 < 1) 
and Figure 4 (f > 1) ; here 2", corresponding to collapse of the mem- 
brane, is greater than 2’, corresponding to rupture of the membrane. 
At any given point in one of these Figures, the sign 2’ is opposite to 
its sign at the corresponding point in Figures 1 and 2. Figures 3 and 
4 show that a cell whose membrane is such that dV /dz < 0 must either 
explode or collapse. 

Finally, a few remarks will be made about the relation between 
the function V (z) and the tension of the cell membrane. If the shape 
of the cell is independent of the volume, conservation of energy dur- 
ing a small change of volume implies that 


0 


P,—P,)dV = 
(P; — Po) ae 


(z—y)dV =a(V)dA=—6a0(V) V4 AV, (28) 


where o(V) is the tension in the cell membrane, in general a function 
of the area and hence of the volume, and 6 is defined as 


Sey 

6=- 
2 A 
a constant if the shape of the cell is invariant. For a spherical cell, 6 
has the value 4(47/3)-7. Equation (23) states that the work done 
by hydrostatic pressure in changing the volume of the cell is equal 


to the work done against the membrane tension in changing the area 
of the membrane. It implies that 


2bK6 
2—y=—— VV a(V). 
Dy 


? 


It seems reasonable to suppose that o(V) is either constant (as in the 
case of the surface tension of a drop of liquid) or else a monotonically 
increasing function of V. If «(V) is the n-th power of V, then, when 
z=y, V is either zero or infinity depending on whether n is greater 
or less than 4; thus, if dV/dz > 0 (n > 4), the cell is collapsed when 
z= y, or every such cell must collapse if it does not burst first, while 
if dV/dz < 0 (n < 4), the membrane is ruptured when z= y. There- 
fore, c= V" is not a possible condition for a durable cell. The quantity 
o«(V) must, in addition to being an increasing function, have a root 
at some value of V, V. = V(y) # 0; furthermore, o must increase 
rapidly enough so that dz/dV > 0, if the cell is to survive. A simple 
possible form for «(V) is C(V — Vo)”, where n is an odd positive in- 
teger and C is a constant. The rapidity with which «(V) increases is 
a measure of the rigidity of the cell membrane, as is the reciprocal of 
the quantity y used above. 
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THE MECHANISM OF THE MIDDLE EAR: 
PART II. THE DRUM 


MARTINUS H. M. ESSER 
ILLINOIS INSTITUTE OF TECHNOLOGY 


The ear drum is considered to be a thin circular membrane with 
radial and circular fibers. whose center is pulled inwards by the handle 
of the hammer. It is shown that such a membrane is equivalent to a 
rigid piston connected by a lever to the handle of the hammer, and sub- 
jected to elastic forces. The stability of the equivalent system is great, 
and the flexibility of the lever is very small. The lever is such that smail 
pressures in the auditory canal are transformed into larger forces on 
the hammer. The leverage ratio increases with the tension of the ten- 
sor tympani and decreases with the number of circular fibers. 


Method. We shall proceed as follows: We first derive the equa- 
tions of equilibrium of the drum, and find the shape of its meridian, 
when subjected to a force ¢ on its center and a uniform pressure p 
on one of its sides. We consider next a small change A¢ of the force and 
Ap of the pressure, and study the corresponding deformation of the 
drum and the variation of its tension. Neglecting all infinitesimals of 
second order, we obtain linear relations between the various incre- 
ments A. From these relations we eliminate all variables relating only 
to internal characteristics of the drum, and so obtain equations re- 
lating the force on the hammer and its motion to the pressure and 
motion of the air in the auditory canal. These equations are: 


Ao =— KnR?Ap — AAX, 
(1) 


They show that the force A¢ on the handle of the hammer de- 
pends linearly on the force 2R?Ap exerted by the air on the drum, 
and on the displacement AX of the hammer, and that the displace- 
ment AX depends linearly on the average displacement AV /nk? of 
the drum, and on the force zR?Ap transmitted. The equations () 
indicate that the drum is equivalent to a system formed by a rigid 
piston, a flexible lever, and a spring which opposes displacements of 
vhe lever from its position of equilibrium. The quantity K is the ratio 
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of the two lever arms, A is the stiffness of the spring, and B meas- 
ures the flexibility of the lever arms. 

Under very particular assumptions, which were originally made 
by H. Helmholtz (1873), we find easily the value of the coefficient 
K, and find also that B = 0. Under more general assumptions how- 
ever, the coefficients K , A, B appear first as complicated expressions, 
and several sections of this paper will deal with the simplifications 
of these expressions. 

Anatomical description of the drum. The ear drum is a thin cir- 
cular membrane (its thickness is 1/10 mm and the ratio of its small- 
est radius to its largest radius is about 4/5). It is tightly stretched 
over a perforation of the temporal bone, and its center is pulled in- 
ward by the handle of the hammer ossicle, which itself is pulled in- 
ward by ligaments and by a muscle, the “tensor tympani”. The in- 
tensity of this pull varies, mainly in animals, with the loudness of the 
sound transmitted. The handle may be considered to be attached to 
the drum at one point because, although in reality the handle is at- 
tached to the drum all along one radius, the pulling force by the han- 
dle seems to be concentrated at the center of the drum. The drum con- 
sists of two layers of tendon-like collagenous fibers and is covered on 
both sides by an epithelial layer. The fibers on the side of the outer ear 
are radial and. the fibers on the side of the middle ear are circular. 
The radial fibers are curved, which indicates that the circular fibers 
are tense. The circular fibers are most numerous towards the pe- 
riphery. 

A-small region of the drum, called “pars flaccida’, lacks fibers 
and is therefore slack. Its existence implies that the fibers do not 
have a strictly radial, circular arrangement and that the air pressure 
in the auditory canal is not always exactly what it would be if the 
entire drum were under tension. The effect of the pars flaccida will 
be neglected in our calculations. For a further description of the drum, 
see, for instance, F. R. Bailey (1936, p. 573). 


Notations. We shall use the following notations: 


= abscissas, measured on the axis of revo- 
lution of the drum, with positive direc- 
tion toward the outer ear, 


r= ordinates of the meridian of the drum, 
or radii of the circular fibers, 


a=a(7r) =angle, at a variable point, of the meri- 
dian with the z axis, 


=s(r) =arc length, 
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R= radius of the circle of attachment of the 
drum to the temporal bone, 
X=a(R) —2x(0) ,S=s(R) —s(0),m=a(0), 


yr =y,(r) = tension in the direction of the radial fi- 
bers, 


Yc = y-(r) = tension in the direction of the circular 
fibers, 


[=I (r) = 2ary,(r) = total tension, along one parallel, of all 
the radial fibers, 


= pull of the hammer on the drum, 


p = excess of air pressure in the outer ear 
over the pressure in the middle ear. 


AAA 


MIDDLE EAR 


R OUTER EAR 


MERIDIONAL 
CROSS SECT- 
ION OF DRUM je- x 


SE eee wali tie aS tes Tet 3 


Equation of equilibrium of the drum. We consider the portion 
of drum limited by a circular fiber of radius r. This portion is in 
equilibrium under the pull ¢ of the hammer, the resultant I’ cos a of 
the radial fibers along the parallel 7 and the force ar’p exerted by the 
air. We have, therefore, 


¢+arp—=—TWI cosa= 2ary, cosa. (2) 


We consider an elementary surface of the drum contained be- 
tween two circular fibers of radius r and r +,dr and between two 
meridional planes forming an angle dé. We obtain two equations of 
equilibrium by writing that this element of the drum is in equilibrium 
first with respect to forces tending to move it in the tangential direc- 
tion and second with respect to forces tending to move it in a direc- 
tion normal to the surface. The two forces acting tangentially on the 
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two arcs of parallel are y,rd# and (yr + dy,) (r + dr) dé = yer a 
d(y,r)d@. The two forces acting on the two arcs of meridian are 
yds. They form an angle 1 — dé, and have, therefore, a resultant 
equal to y.ds dd. This resultant is directed along the radius. Its pro- 
jection on the tangential plane is y. ds dé sin a = y, dr dé. Since this 
projection equals the tangential resultant d(y,r)dé of the forces on 
_the ares of parallels, we obtain 


_ ad(ry) 1 @P 
Se dp ae dr- 


Ye (3) 
The condition of equilibrium with respect to the forces normal to 
the surface leads towards the known equation 


Vr Ve 
| ee Smears ? 

Pr Pc 
where p, and p- are the two principal radii of curvature of the sur- 
face. These radii are positive if the corresponding centers of curva- 
ture are in the region containing the fluid under pressure p , and nega- 
tive in the contrary case. We have 

ds dr dr 
pr _— = 


da sinada dcosa_ 


Two normals to the membrane at two points of a same parallel inter- 
sect on the axis of revolution. Therefore p, = 7/cosa. We thus ob- 
tain the equation 


d cos a COS a 
c 
dr ? 


The three equations (2), (3), and (4) are not independent, and equa- 
tion (4) is a consequence of equations (2) and (3). 

Equation of the meridian. We have thus far two independent 
equations, namely, equations (2) and (3), to determine the three un- 
known functions r(x), y,(r), yc(7). A third equation is required. 
This equation depends on the relative distribution and tension of the - 
circular fibers, and must therefore be obtained by considerations other 
than purely static. 

If the shape of the drum is given by an equation r = r(x), then 
we obtain successively a(r) by tan a= dr/dx, y,(r) by equation (2), 
and y.(’) by equation (3). 

Conversely, if the function y,.(”) is given, then we obtain yr(r) 
from equation (3), which we write in the form 


P=Yr (4) 
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.T r 
Ty, = —= { Vo(r) ar + a, (5) 
JU 0 


a being a constant of integration. The function a(r) is then obtained 
from equation (2) and the function x(7) is obtained from 


T r 


veers (¢ + apr?)dr | 
x(r) = i =| : (6) 
tan a Vid +i pr) 2 


0 0 


If y, is given as a function of s instead of as a function of 7, and 
if p = 0, we use the following expression: 


ios easy d d 


onde 2h sin a ds cos a 
pS 1 ad ade 1 da__@ dtana 
2a sin acos?a ds 2ncos?ads 21 ds — 
Therefore, 


Qn s 
ace a, [ vel) ds + tan ay. (7) 


We thus obtain tan a = dr/dz as a function of s. The functions x(s) 
and r(s) are then obtained from 


s s 


ids tan a 
i ———————_, r= —__—_—— ds. 
V1 + tan? a \/1 + tan? c 


0 0 


Case where y. is constant and p is zero. We shall apply to this 
particular case the results of the preceding section. Defining two con- 


stants b and c by 
b=Y¥e5 c= ¢/(2ab), 


the equations (5), (2), (7), and (6) give respectively: 


if be 
Were = RENE ar OT ea Qo = CS, 


bedr a+ br |" 
= | = ¢ arg cosh < 
a (oro 0 ok all 


0 
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Thus the meridian appears to be a hyperbolic cosine curve. For such 
curves we have the following known relations, d being a suitable con- 
stant: 


a Cri O, ee, 
r +-=c cosh —— = V2 1+ e= ’ 
b G COs a 


dr +d 8s 
een === tamer (8) 
dx G Cc 
a a \2 od. 1+ sina 
log BRN (elle — log A =—— = log — ; 
b b Cc COS a 


The condition that these equations have a meaning for r = 0 requires 
that a 2 bec = ¢/2a. 


The case where we take as the third equation y.(r) = y,(7) in- 
stead of y.(r) = b is of interest because the condition y, = y, is char- 
acteristic of amorphous surfaces or surfaces with random disposition 
of fibers (such as rubber membranes) and also of surfaces of mini- 
mal area. In that case, equation (3) gives y. = y, + dy,/dr. There- 
fore, the equation y, = y, implies that y,, and therefore also y., is a 
constant. It follows that equations (8) are again valid, and that 
a=0. Since for a= 0 equations have a meaning only for r 2 ¢, it 
is necessary to suppose that the attachment of the hammer to the 
drum is circular and has a radius at least equal to c. 


The hypotheses of Helmholtz. H. Helmholtz (1873, appendix) 
makes the following three hypotheses: 


1. The shape of the drum is the one it would have if there were no 
circular fibers but if there were instead an additional pressure q 
in the middle ear. 

2. The radial fibers are inextensible, which means that their length 
is independent of their tension. 

3. The circular fibers are extensible, which implies that their ten- 
sion changes only when their length changes. 


The interest of these hypotheses will be seen in the following 
section. In the present section we shall determine the functions a(r), 
x(r), yv,-(v) and y-(7) resulting from the first hypothesis. 


From assumption (1) we obtain a(r) and x(r) by setting ve = 0 
in equation (5) and p = —q in equations (2) and (6): 
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ct Pa eae 
cos a(7) =A 7s COS ao (1 — ar?q/¢) , 


i (9) 


x(r) =| erent —dr. 
V (22a)? — (6 —a qr’)? 


0 


This is an elliptic integral. H. Helmholtz (1873, pp. 67-69) studied 
this integral and gives a figure of the meridian for a, = 40°. 
Actually, however, y, # 0 and p= 0. Helmholtz’s first hypothe- 
sis amounts to the assumption that equations (9) hold. We obtain 
yr(r) and y.(7) from equations (2) and (3): 


eee eee CT 
COSa CoOSa 1—arg/d 


SS A RAR qr ou 


Ce ae Oo (1—27°q/s)? 


Ve 


The last equation agrees with the observation that the density of the 
circular fibers increases towards the periphery of the drum. 


Transmission of forces under the Helmholtz hypotheses. The 
Helmholtz hypotheses lead towards some simple conclusions, because 
they enable us to study different conditions of equilibrium without 
considering deformations of the drum. We shall show that the posi- 
tion of equilibrium of the drum does not depend on the pressure in 
the auditory canal as long as there is no displacement of the hammer. 

Let us consider a pressure Ap in the auditory canal. Leaving the 
‘function a(7) unchanged, we increase all stresses and the pull on.the 
hammer by amounts equal to — Ap/q times their values when there is 
an excess of pressure g in the middle ear, but no circular tension. Then 
the equations (2), (3) and (4) of equilibrium of the drum, being 
linear and homogeneous in y,, ye, ¢ and p, remain satisfied for the 
new state of the drum, and the equations of elasticity of the fibers 
are also satisfied, for we have no change in the length of the radial 
fibers, and no change in length and tension of the circular fibers. We 
have, therefore, a new state of equilibrium. The first Helmholtz hy- 
pothesis was made in order to have the increments of the different 
forces satisfy the linear equations (2), (3) and (4), the second hy- 
pothesis in order to avoid a stretching of the radial fibers due to the 
- increase of their tension, and the third hypothesis was necessary in 
order to have the values of the tensions determined uniquely. 
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In particular, the increment 4¢ of the pull ¢ in the new state is 


(peat Ly rs ? nah? Ap. (11) 
q nk? q 


Thus to a force 2R?Ap exerted by the air on the drum corresponds 
a force —A¢ exerted on the hammer, which is ¢/(zk?q) times larger. 
Therefore, for transmission of forces without motion, the drum is 
equivalent to a piston connected to the hammer by a lever increasing 


forces in the ratio 


ull 
¢ ~*~ forceot pul (12) 


ak?q force of pressure 


This lever is rigid [B = 0 in equation (1)]-because for an immobi- 
lized hammer there is no motion of the drum whatever the force act- 
ing on it. The ratio K can be made as large as we want by making 
q, and therefore the circular tension y, [cf. equation (10)], small 
enough. The ratio K can in particular have the value 3 or 4, which 
was found to be most suitable for the ear (Esser, 1947). The ratio 
K depends on #. Therefore the efficiency of the ear can be modified 
by changing the tension of the tensor tempani. 


Displacement of the drum; discussion. of the Fane Pi The 
preceding sections gave the essential properties of a non-moving drum. 
In the following we shall study displacements of the drum from its 
position of rest. No assumption will be made about the distribution 
of the circular fibers. In other words, the function y,(7) will remain 
indeterminate. We shall, however, suppose that the circular stress 
y-(7) does not change during the motions of the drum. This assump- 
tion is made for mathematical simplification, but is not justified by 
physical or anatomical considerations. It would be justified 


1) if the elongation of the circular fibers were small in comparison 
to the elongation of the radial fibers (movements parallel to the 
x axis), 

2) or if the circular fibers were very stretchable, in the sense that 
small changes in tension were to induce large elongations. 


Our hypothesis is to some extent similar to Helmholtz’s hypoth- 
eses (2) and (3) in the sense that it supposes the circular fibers much 
more extensible than the radial fibers, but with the emphasis on the 
extensibility of the circular fibers rather than on the inextensibility 
of the radial fibers. As far as we know, no such difference of elas- 
ticity of the two sets of fibers has been observed in the living ear. 

We shall use the notations y(r), x(r), a(r), 6, p = 0 , ete., to 
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denote functions and quantities relative to the rest position and will 
designate the same functions or quantities by y(r) + Ay(r), x(r) 
+ Ax(r),a(r) + Aa(r), ¢ + Ad, Ap, ete., at an arbitrary moment. 
We suppose that all increments A are small enough to allow us to neg- 
lect terms of the second order and therefore we can apply all the rules 
of differential calculus to the increments A. 
As we have p = 0 and Ay. (r) = 0, we obtain from equations (2) 
and (5): 
Tr csa=¢, (13) 


A I = 2nAa = constant. (14) 
Thus AI is independent of 7. 


Displacement of the hammer. In this section, we shall calculate 
the displacement — AX of the hammer in terms of 4d, Ap and ATI. 


R 
- We have X = | dr/tan a, therefore, 
0 


R 


1 
ax=[ dr. (15) 
tan oa 
0 
We have also 

1 A A cosa 
SS ens (16) 

tana sin? a sin? a 


and, from equation (2), we obtain 
otarp cosa 
r 


We introduce the following notations, some of which will be used now 
and some of which will be used in later sections: 


cos”? a 
Acosa=A (Ad + ar?Ap) LP EIE: heres (17) 


R R R 
iy 
COS a cos? a ) cos® a 
ifise dr, MS = dr. = og dr 
sin? a sin? a sin? a 
s 
0 


Rk 
Me 2 
fe Bd) Eee af hea (18) 
R2"} sin’ a ’ R2 | sin?a 


84 THE MECHANISM OF THE MIDDLE EAR 


Expressions (18), together with equations (15), (16), (17), give 
R? ADT 

Omid oy ar (19) 

p oo) ¢ 


Calculation of arc length and volume. We calculate now the in- 
crease AS of the length of the radial fibers in terms of 4¢, Ap, Al’. 


R 
We have S = i] dr/sin a, therefore, 
0 


1 
AS = [ dr, 
sin a 


(20) 
1 ..cos 4 A cos a 


sin a sin? a 


Thus the integrand in equation (20) differs from the integrand in 
equation (15) by a factor cos a. Therefore, AS is obtained by replac- 
ing in the second member of equation (19) the integrals L, P, M, 
by M,Q,N respectively. We thus obtain 


AS =— M + ——Q—-—N. (21) 
i) ? co) 

We now calculate the increase AV of the volume of air contained 
in the auditory canal. The quantity A4V/(zR?) measures the average 
displacement of the drum or the displacement of a rigid piston equi- 
valent to the drum. We have 


x R 


v=a | r? (a2) toxn [ Bay, oe 
tan o 


0 0 


R 


1 
av=s { v4 aie 
tan o 


0 


The integrand in this equation differs from the integrand in equation 
(15) by a factor 7?. Therefore, AV/z is obtained by replacing in the 
second member of equation (19) the integrals L, P, M by R°P, R°T, 
RQ respectively. We thus obtain 
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AV Ad nk? Ap Ar 
Bees teat T ——— Q. (22) 


Elimination of AI. We shall find the relation existing between 
AS , AI and the coefficient k of elasticity of the radial fibers. This 
relation will enable us to eliminate AI’ from equations (19) and (22). 

The coefficient & of elasticity of a fiber is defined by Ay/ y 
= kAs/s , where y indicates the tension at rest of an element of fiber 
of length s and y + Ay indicates the tension of the same element when 
- it is stretched to the length s + As. 

The increase AS of the length of a radial fiber for a given value 
of AI is calculated as follows: 


1 Ay, AT ds 
AS = -—— 3 = —— — 
ky, k r 
Ss xi 
AT A 
_ AD cosa Ar B 
k ¢ 
x 
AS=— Ar. (23) 
ke 


This is the desired relation between AS and AI’. Elimination of 4S 
from equations (21) and (23) gives 


A R?A 
= (n+2)=2u+o. (24) 
d k 

Substitution into equation (19) gives 


ax (w=) =*[21+Nz—ar| 
k dLk 


= ak (=P NP—mQ| 
Solving the preceding equation for 4¢, we obtain 
Ao =— KaR?Ap + AAX, (25) 


where aay 
__ (X/k)P +NP—MQ> Re ea aE A (26) 


 (X/k)L + NL— M2’ ~— (X/k) L+ NL—™? 
This proves the validity of the first equation (1). 
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Elimination of 4¢/¢ and AI'/¢ from equations (19), (24) and 
(22) gives, in the form of a determinant, 


Cem — AX + Pak? Ap/¢ 
M N+X/k Qrk?Ap/d| =0. 
PO — AV /aR? + TnrR?Ap/o> 
Expanding this determinant, we obtain 
1 AV 
: AX = ———— BaR?Ap, (27) 
Kank one 
where K has the value given by equation (26) and where 
CoM ae 
ENE) ee eee) (28) 
(X/k)P+NP—MQ ¢ Pig 7 


By establishing equation (27), we prove also the validity of the 
second equation (1). An interpretation of these equations (1) in 
terms of an equivalent mechanical system has been given in the first 
section. It should be remembered from our definitions that AX or 
AV are positive when the hammer or the drum move toward the mid- 
dle ear, which is in the negative x direction. 

Introduction of multiple integrals. In the coefficients K, A, B 
in equations (25) and (27) appear the four determinants: 


Dearne 

N M LM for 

Q P le / |p 7 [ana D=| % N Ql - 
PQ TI (39) 


We shall transform the three first determinants into double integrals 
and the determinant D into a triple integral. These multiple integrals 
will be used in later sections to find the order of magnitude of the 
determinants for small values of a(R) — a(0), to calculate their ap- 
proximate values and to find the range of possible values of K . 

We define a function u(7) by 


r 


COS a 
u(7’) ={ - dr. (30) 
sin? a 


0 


This function u(7) is monotonically increasing, and therefore the in- 
verse function r(u) is well defined and monotonic. We use also the 
notation 


c(u) =cos alr(u)]. (31) 
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Then the integrals (18) become 


L oL L 
L= fw, m= | En aus N= [ ct(u)dw, 
0 0 0 


L L L se 
nares r2 d Q Suk v2 3 d d pe A 4 d: 
S (uw) du, Saye (u)c(u)du, apr ee (u)du. 
We have 
[ fte@—ew)  , 
0 (0 (33) 


- [e(u) 7? (v) — c(v) r?(u)] dudv = 2R? (NP — MQ), 
because, multiplying out the integrand, we find 


L L L L 


f qa { e@)ao— J e(uyre (ua f e(eyaw 


L L 


= Recs J ecoyr2(a) dv + fra f e@yae 


= N (PR?) — (QR?)M — M (QR?) gle (PR?) N = 2k? (NP — MQ). 


By a similar argument, we obtain the formulae 


2(LN—Mty= | [ ec e(us lsdadee (34) 


2R*(LT — P2) = [ f te@ —K eK dudo, (35). 
Rp= | { f Letu) — eo) ILe(v)r#(w) — e(w) 7? (0) 
0-0 0 (36) 

. [r?(u) — 7? (w) ]du dv dw. 


Thus we have expressed in equations (33), (34), (385) and (36) 
the determinants (29) as multiple integrals. The integrands in equa- 
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tions (33), (34) and (35) are positive; therefore the first three de- 
terminants (29) are positive. 

Magnitudes for small a(r) —«(0). As in the case of Helmholtz’s 
hypotheses, it will be found that a(r) — a(0) must be small in order 
’ to obtain a large value of K. When a(r) —a(0) is small, the determi- 
nants (29) appear as differences of approximately equal terms, and 
are therefore inconvenient to use. However, the multiple integrals 
are of such a form that their magnitude can be estimated. 

The integrands in equations (33), (34), and (385) have the same 
(positive) sign over the whole domain of integration. Therefore, the 
corresponding integrals have the same order of magnitude as their 
integrands. For small values of a(r) — a(0), the difference c(u) — 
c(v) is an infinitesimal of first order, while c(u)7r?(v) — ce(v)r?(u) 
and r?(u) — r?(v) are not small. Consequently, LT — P? is not small, 
NP — MQ is an infinitesimal of first order, and LN — M? is an in- 
finitesimal of second order in a(R) — a(0). 

The determinant D is an infinitesimal of at least first order in 
a(R) — a(0), but as the integrand in equation (36) does not remain 
positive, D may be an infinitesimal of higher order. We shall show 
that it is at least of second-order. 

We have 


fifefteco —e(v) ][72(w) —72(v)] 


- [7?(u) — 7? (w)] dudvdw=0, GP 
because the integrand is an odd function of «— v and the domain of 
integration is symmetric in u and v. We multiply the two members 
of the preceding equation by an arbitrary constant A and subtract 
from equation (36). We thus find that equation (36) remains valid — 
when we replace the integrand by 


[e(u) —e(v)][{e(v) — A} r2(w) 
LCM s— A yrete) |r (a) ete 


If we now take for A an average value of c(v), it becomes evi- 
dent that the integrand is an infinitesimal of at least second order in 
a(R) —a(0). 

We shall now find the range of possible values of the ratio K , 
defined by equation (26). The denominator of K is positive for all 
values of the coefficient of elasticity k. As k varies from zero to in- 
finity, K varies continuously and monotonically from the value P/L to 
the value -(NP — MQ)/(NL — M?). The value P/L approximately 


/ 
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equals 1/3, as is seen by taking w proportional to r in formulas (32). 

The value (VP — MQ) /(NL — M?) can be made as large as we want 
by taking a(R) — a(0) small enough. Therefore, the range of pos- 

sible values of K is 1/3 to infinity; the smaller values of K being ob- 
tained for small values of k , that is, for very stretchable radial fibers; 
and the large values of K being obtained for large k and small a(R) 
— a(0), that is, for unstretchable and tense radial fibers. 


Computation of the integrals.. To compute the integrals, we shall 
suppose that cos o(7) varies proportionally to the square of 7: 


cos a(7) =cos a(0) — [cos a(0) — cos a(R) ] (7/R)?. (38) 


This relation holds when the first hypothesis of Helmholtz is satisfied 
[cf. equation (9)]. Using the notations of equation (31), we have 


c(u) — e(v) = [cos a(0) — cos a(R) ]R? [7?(v) — 7? (u)], 


c(u)r*(v) — e(v)r? (vu) = cos a(0) [77 (v) —7?(u)], Se 
and the integrals (33) and (34) become \ 
2R? (NP — MQ) —cosa(0) [cos a(0) — cos a(R) ] 
-R? ff [7r?(v) —7?(u) ]? du dv, 
2(LN — M?) = [cos a(0) — cos a(R) ]? 
-R-+* f§ f[r?(v) —7?(u)]? dudv. 
Using also equation (35), we obtain 
NP — MQ = cos a(0) [cos a(0) — cos a(R)](LT — P?), ns 


LN — M?= [cos a(0) — cos a(R) ]?(LT — P?). 


Furthermore, to calculate the integrals L, N, P of equations (18) 
and LT — P? of equation (35), we shall neglect the curvature of the 
radial fibers; in other words, we shall assume that a(7) is constant: 


cos a cos? a lecosa 


sinta ” sinta ” 3 sin? a 
(41) 


hy cos? a 7 4 cos?a " 
= — s?)?drds=———_ BR 
LT — a 2R* sin’ a ‘f f Seas 45 sin’ a 


From equations (40), then, we obtain 
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4 cos® 
NP —MQ—— ——~ [cos «(0) — cos a(R) ] R?, 
45 sin*® a 
(42) 
4 cos? 
LN — M?=———~ [cos a(0) — cos a(R) ]? R?. 
45 sin® a 


We still have to calculate the integral D of equation (36). Using re- 
lation (39), we obtain 


R*D = cos a(0) f f fLe(u) —e(v) [7 (w) — 7? (v)] 
-[r?(u) —7r?(w)] du dv dw. 


Therefore, by equation (37), we have 
D=0- (43) 


Remark. The computation of the integrals (42) and (43) was 
based on the assumption that a(7) varies as the square of r. The val- 
ues obtained, however, differ little from the values obtained with other 
hypotheses. For instance, if we had supposed that the radial fibers 
are arcs of a circle, which is to say that cos a(7) varies as the first 
power of 7, we would again have obtained the equations (42), except 
that the coefficient 4/45 would be replaced by 1/12. The value of D 
would not be zero, but 

D=—— 2 Fos a(0) (R)}? Rs 
=—— — [cos =— CO. gl Oe 
2160 sina ee 

Computation of the coefficients K, A,B. We have approximate- 

ly R=X tana, and 


cos a(0) — cos a(R) =sin a[a(R) —a(0)]. 


Using these relations and equations (41), (42), and (43), we ob- 
tain by substitution in equations (26) and (28): 


15+ 4 cotala(k) —a(0)] & 
~ 45 + 4[a(R) —a(0)]2k 

_ 45sin?a+ 45cos*ak sinad 
~ 45 + 4[a(R) —a(0)]2k cosaR’ 


? 


(44) 


pe 4 . \,€0s 0K 
15 + 4 cot a[a(R) —a(0)] k sintad 
Moreover, it may be convenient to express a(R) — a(0) in these for- 


mulae in terms of the tensions in the drum, as we shall do now. We 
have from equations (7) and (13), 
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8 


a(R) —a(0) 2 
——__——- = tan a(R) — tan a(0) =——— | ».(s) ds. 
COS” a I’ cosa 
Therefore, 
So ye(s) ds 
a(R) —a(0) = 22cos a —————_, (45) 
Wg 


where f* y-(s)ds is the sum of the tensions of all the circular fibers 


and I is the sum of the tensions of all radial fibers. 

If we suppose the radial fibers unstretchable (k = o) and that 
sin a= 1, and if we use equations (45) and R= X tana, we can sim- 
plify the equations (44) into 

; oho ieee Pps B=0. (46) 
2nf* y-(s)ds 4X 

Conclusion. The dynamics of the ear drum lead to the equations 
(1), where the coefficients K , A, and B are given by equations (44) 
if we suppose the radial fibers stretchable, and by equations (46) if 
we suppose the radial fibers unstretchable. An interpretation of these 
equations has been given earlier. 
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A MATHEMATICAL DESCRIPTION OF METABOLIZING 
SYSTEMS: II 


HERMAN BRANSON 
HOWARD UNIVERSITY 


Some of the laboratory procedures available for determining the 
functions in the integral equation established in part I are discussed. 
The tracer or tagged molecule technique is shown to be especially prom- 
ising including the use of “double tracer” molecules. Conversely, the in- 
tegral equation may be a convenient device for correlating and integrat- 
ing some of the work now being done with tracer molecules in biologi- 
cal systems. 


In part I of this series of papers (Branson, 1946), the author es- 
tablished the integral equation 


M(t) =M(0) F(t) + [ro F(t—6) do (1) 


to describe the behavior of a metabolite, M(t), in a system with rate, 
R(t), and metabolizing function, F(t). The usefulness of this equa- 
tion will depend upon the experimental techniques available for de- 
termining any two of the functions in order that the third may be 
derived by solving the equation. 

Although there are many laboratory techniques which will deter- 
mine M(t) as a function of ¢ and control R(t) so that F(t) can be 
determined for certain specific conditions, the author is especially in- 
terested in determining these functions by the use of tagged or tracer 
atoms (radioactive or stable isotopes of ordinary chemical constitu- 
ents of the metabolite). In all experiments, however, it is usually 
possible to plot M(t) vs. ¢ and determine by the usual mathematical 
techniques (Worthing and Geffner, 1943) a continuous function which 
will approximate the data within the experimental error. 

The Determination of the Functions in the Integral Equation: In 
the typical experiment on metabolism, we are interested in the fate 
of some metabolite which we shall designate as M(t). If our experi- 
ments are on mature, normal animals in nutritional equilibrium, we 
often find that the amount of the metabolite varies little if at all. The 
introduction of a tagged quantity of that metabolite, M(t), reveals 
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that the constancy is only apparent. The equilibrium is dynamic, as 
has been so ably demonstrated by R. Schoenheimer and collaborators 
(Schoenheimer, 1946) with the stable heavy isotopes of hydrogen 
(H?) and nitrogen (N°). Hence when we assay for the metabolite, 
we are on the flat portion of a curve whose origin cannot be deter- 
mined by examining M(t) at the time of our experiment. According 
to equation (1), we have for this situation 


M(0) [1—F(#)] =f R0) F(t—6) do, (2) 


with R and F still to be determined. We need another independent 
relation for determining one of them. 

The tracer atom technique gives us the required condition, for we 
need only inject a sample of tagged metabolite, M*(0), and follow its 
course in the same system or systems sufficiently similar. These data 
will give 

M*(t) = M*(0) F(é). (3) 


Since no additional tagged metabolite can enter, R* = 0. Thus, by 
means of equations (2) ‘and (38), we can determine the functions & 
and F. 

The preceding formulation of equation (3) assumes that al- 
though the rate at which the tagged substance enters is different from 
that of the normal, the metabolizing functions are identical for both. 
If that were not true, then the system would be discriminating be- 
tween the tagged and untagged molecules of otherwise identical chemi- 
cal constitution. Within the experimental error of this work, the evi- 
dence is convincing that the biological system cannot distinguish the 
small differences in mass between the chemically similar units. There 
may exist slight differences in rates of diffusion and other physical 
phenomena. 

If we are dealing with a metabolite or substance which is not in 
equilibrium, the procedure is the same except that then we must de- 
termine M(t) and use equations (1) and (3) for R and F. 

Treatment of some Data: As examples we shall consider some 
experiments using radioactively tagged molecules of J. G. Hamilton 
and M. H. Soley (1940), some experiments and theoretical develop- 
ments by D. B. Zilversmit and collaborators (1943), and we shall al- 
lude to some experimental results of R. Schoenheimer and his collab- 
orators (Schoenheimer, 1946). 


J. G. Hamilton and M. H. Soley (1940) fed radioactive iodine 
(I+) to a group of normal human subjects and observed by means 
of an external Geiger-Miiller counter the emanations from the thyroid. 


a‘ 
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From the many types of empirical curves which may be fitted to their 
data, the writer selects 


M(t) =C(1—e*) oF! 


with C = 0.035 M*(0), where M*(0) is the amount of tagged iodine 
fed. This choice of M(t) describes the observed behavior. satisfac- 
torily and the expression is sufficiently tractable mathematically for 
straightforward integration in equation (1). From their data we see 
that the maximum of the curve is reached in about one day, hence 
a © 4.5 days?. After 30 days, over 80 % of the radioactive iodine re- 
mained in the thyroid, hence 8 ~ 0.006 day. Inasmuch as there was 
no radioactive iodine originally present, equation (1) now becomes: 


Gil — e-#').e#' = ez F(t—6) dé. 


Although we cannot write an equation such as equation (3) on 
the basis of the data, the fact that the radioactive iodine slowly leaves 
the thyroid limits the possibilities for F(t). Upon integration and by 
use of this qualitative information on F(t), we have 


F(t) = et, 
R(t) =a eo", 


Thus we may conclude on the basis cf our formulation that the 
metabolizing function describing the history of iodine in the thyroid 
is a slowly decreasing function of time, while the rate function de- 
creases rapidly with time. Although the parameters « and f are em- 
pirical and are not related to physiological processes, nevertheless, 
their values may be important clues in detecting malfunctioning of 
the thyroid. 

The example given by D. B. Zilversmit and his collaborators (Zil- 
versmit, Enteman and Fishler, 1943; Zilversmit, Enteman, Fishler 
and Chaikoff, 1943) is an example of the general problem of the con- 
version of A into B where A is called the precursor of B. A may be 
produced by complex reactions and B may be lost through others. We 
shall assume only that for each unit of A which disappears a unit of 
B appears. If A and B are tagged substances, we shall have none of 
either present initially, hence 


(4) 


A*(t) = {Rw F(t—6) dé, 


B*(t) = [R.©) F,(¢—6) do. 


96 METABOLIZING SYSTEMS 


If we assume that the transformation of A into B is an ir- 
reversible, first order reaction, chemical kinetics states that R,(0) 
= CA*(6), whence 


B*(t) = 65| v7 {Rw F(6@—¢) ap | F.(t—9) dé. 


The resulting expression for B*(t) is complicated as one would 
expect in the absence of additional simplifying assumptions. There 
are three functions to be determined, R, F and F,. We shall need 
one additional relation in order to determine the functions uniquely. 
One approach would be to follow the behavior of A* and B* in a simi- 
lar system where it is injected. For such a system A*(t) = A*(0) F(t). 

The systems treated by D. B. Zilversmit and collaborators (Zil- 
versmit, Enteman and Fishler, 1943; Zilversmit, Enteman, Fishler 
and Chaikoff, 1943) are simpler. Their general system (A > B) is 
described in our integral equation formulation by the equations 


A*(t) = [ R) F(t— 6) do, 
B(0) =B(0) F.(t) + J Ra(8) F(t—9) a0, 
Be (6) 6 {'ar(e) F,(¢— 0) dé, 


which state that the amount of B present is constant, and R = CA 
which is constant if the amount of A is constant. Thus from follow- 
ing the courses of A*, B*, and determining B in one system and either 
A in the same system, or A* or B* in a similar system, we shall have 
enough relations to determine the R’s and F’s. If we are interested 
only in the & and F associated with B, we need know only A*(¢t), 
B*(t), and B(0Q). 

Their experiments on the turnover rate of phospholipids in the 
plasma of the dog with radioactive phosphorus follow the conditions 
for the system described by equation (9) of part I for the ordinary 
phosphorus and B*(t) = B*(0) F(t) for the radioactive. Since R 
is constant, F(t) = exp (—R/M)t. 

The curves describing the atom percent deuterium uptake of 
cholestorol in mice (Schoenheimer, 1946), as well as other experimen- 
tal results with stable isotopes (Schoenheimer, 1946), are expressed 
with acceptable accuracy by M(t) = C(1 — e**) e8', witha > > B, 
which gives expressions similar to equations (4) for the metabolizing 
function and the rate. 

The application of double tracer techniques, e.g. using radioactive 
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carbon (C**) and “heavy” carbon (C**), H? and H?, etc., to label the 
molecules would be particularly advantageous in determining the 


functions in equation (1). We could determine simultaneously A*(t) 
in 


A*(t) = { R@) Fa—4) dé 


and A*(t), where the bar indicates the second type of tagged mole- 
cule, in A*(t) = A*(0) F(t). Thus by using two differently tagged, 
but chemically similar, molecules in a single system, we would have 
the two equations necessary to determine both R and F. 

As a specific example of the double tracer technique let us con- 
sider an experiment on the metabolism of methionine in the rat or 
mouse. We can feed homocystine labelled with radioactive sulfur 
(S**) and methionine labelled with a stable rare isotope of sulfur 
(S* or S**). It may be necessary to feed some choline or betaine in 
order that the animals may synthesize the methionine from the homo- 
cystine. From a series of animals, we may isolate methionine from 
the body tissue. Radioactive assays will give A*(t), the radioactive 
methionine synthesized from the homocystine. Determinations with 
the mass spectrometer will give A*(t), the amount of methionine re- 
maining from the labelled amount fed. There are severe experimental 
difficulties in the determination of sulfur; nevertheless, the proce- 
dure seems practicable and such experiments are being planned in 
our laboratory. 

The preceding discussion reveals how the tracer element tech- 
nique may be used in determining the needed functions in equation 
(1). It may be, however, that the chief significance of equation (1) 
will be its value in integrating and correlating some of the work now 
being done with biological systems using the tracer element and other 
techniques. 

The author wishes to express his appreciation to Research Cor- 
poration of New York for a grant to support co-operative research 
in chemistry and physics with stable isotopes and radioactive isotopes 
which has motivated this work. He is indebted to Dr. M. F. Morales 
for a careful criticism of these papers and several fruitful suggestions. 

The publication of this paper has been assisted by the Publication 
Fund in the Natural Sciences of the Graduate School, Howard Uni- 
versity. 
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A MATRIX CALCULUS FOR NEURAL NETS: II 


H. D. LANDAHL 


THE UNIVERSITY OF CHICAGO 


_ ,_ Im a previous paper a method was given by which the efferent ac- 
tivity of an idealized neural net could be calculated from a given affer- 
ent pattern. Those results are extended in the present paper. Conditions 
are given under which nets may be considered equivalent. Rules are 
given for the reduction or extension of a net to an equivalent net. A 
procedure is given for constructing a net which has the property of con- 
verting each of a given set of afferent activity patterns into its corre- 
sponding prescribed efferent activity pattern. 


In a previous publication (Landahl and Runge, 1946), we de- 
fined a structure matrix F = |f;| related to a neural net of N neurons. 
Rach element /;, of F is a measure of the potential action of a neuron 
N; on a neuron N;. The matrix F contains four non-zero sections Fx , 
F,, F; and Fz which determine respectively the relationships affer- 
ent-efferent, afferent-internal, internal-internal and internal-effer- 
ent. We also defined a normalizing operator G and the activity vector 
a(t) =r(t) + i(t) + e(t) which specifies the activity of all neurons 
of the net at time t. The (1 X p) vector r(t) describes the activity 
of the p receptor neurons at time t, the (1 X ) vector i(t) describes 
the activity of the internal neurons, while the (1 X «) vector e(t) 
describes the activity of the efferent neurons at time ¢. An equation 
was given by which any a(t) could be determined from a knowledge 
of the sequence of r’s and a set of initial conditions [loc. cit. equation 
(10) ]. 

We introduce some definitions. Suppose that two nets, with N’ 
and N” (> N’) neurons, each containing the same set of afferents and 
efferents, have the property that if ro(t) = 0 for all t < 0, then for 
any sequence Ro(r) = [ro(0), ro(1), ---] there results a sequence 
E.(e) of é,’s, the same for both nets. Then N’ is a reduced equivalent 
net of N”, and N” is an expanded equivalent net of N’. If the above 
holds, not for every sequence, but only for a given set of sequences 
of r’s, then we may refer to N” as an expanded adequate net of N’ 
over the set of pairs of sequences Reo, Ec. 

Consider two matrices F and F’ of orders (p + « + e) and 
(p + 1 + «) which have no internal neurons in common. The matrix 
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given by their direct sum represents the sum of the two nets if there 
are no neurons in common. On rearranging rows and columns, this 
matrix can be transformed into a square matrix of order p + p’ + 
+ ;'+e+e'. If the two nets have the same set of afferents and effer- 
ents, this matrix can be reduced to a square matrix of order p + 1 
+.’ + « by adding the two rows which represent the same neuron 
and adding the corresponding columns, then deleting the duplicate 
row and column. The resulting matrix shall be referred to as the ex- 
tended sum of F' and F’. It may be written as: 


0 EE; he 0 F'p F's 
0° Fy, 2: 0 es eo a 0 Bs (1) 
eee On 020 ae 

ee aint ae 


where F”y = Fy + F’y. 

If F and F" represent equivalent nets, then their extended sum 
represents an extended adequate net of either primary net. But if 
Fy, = F’y = 0, then the extended sum represents an equivalent net 
if for each BA eolumn of the final matrix it is true that not both 
columns of the submatrices. F’, and F’; contain elements such that the 
sum of any group in the column lies between zero and one. If matrices 
F and #” contain neither finite negative elements nor proper fractions 
in the sections Fy and F,,, or if the sections contain only integers, 
then their extended sum represents an extended equivalent net. If F 
and F” represent equivalent nets, if the simplifying rules have been 
applied to eliminate non-functional elements, and if f,. denotes an 
element of F'; , then f,. and f’,. are both greater than or equal to one, 
both equal to zero, both arbitrarily large negative elements, both 
proper fractions, or are both finite negative elements. If Fy and F’y 
contain no proper fractions or finite negative elements, then Fy = F’y 
and Fy may be replaced by Fx in equation (1). 

We may now write the following efferent separation theorem 
(I): Any net represented by the square matrix F’ of orderp +4 +e 
can be replaced by an expanded equivalent net, the duplicated net, 
represented by a square matrix F" of order p + eu’ + e which is the 
extended sum of F’x, F',,--- F’., --- Fe where F’. is F’ — F’x with 
all of the efferent columns except e replaced by columns of zeros. 

From an inspection of F” it becomes clear that a'(t) and a” (t) 
differ only in that the internal vector c” (¢) is a vector (1 X 1 e) which 
results from extending the (1 X 1) vector u(t) by repeating it «times. 

The following rules may be useful in simplifying computations 
and enabling reductions to be made according to rules given ice 
quently: 
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1. If any column in the structure matrix F contains only one 
non-zero element, this element may be replaced by one if it is greater 
than one, zero if less than one. 

2. If any column in F contains no negative elements except ar- 
bitrarily large ones, then any element greater than one can be re- 
placed by one. 

3. Any non-functional element f;, may be deleted from F. Such 
an element is less than one, and satisfies the condition (S; + #4)G = 
Si G where S;, is the sum of any group of elements except f;, in column 
k of F. 


We next give some rules which result in reduced equivalent nets: 


1. If any row of the structure matrix F representing an inter- 
nal neuron contains only zero elements, that row, the corresponding 
column in F and the corresponding space in the activity vector a may 
be deleted. 

2. If any column of F representing an internal neuron contains 
only zero elements, or if the sum of the positive elements of that col- 
umn is less than one, that column, its corresponding row in F and 
the corresponding space in a, may be deleted. 

3. If any two columns in F representing internal neurons are 
identical, one of these columns and its corresponding row in F' and 
the corresponding space in a may be deleted, provided that the row 
corresponding to the column not omitted is replaced by the matrix 
sum of these two rows. 

4. Suppose that in the matrix F, two rows «a and f are identical 
and that each represents an internal neuron. Consider the following 
set of conditions: 

a) The rows of a and f contain no proper fractions or finite 
negatives. 

b) Neither of the corresponding columns a and f contains finite 
negative elements. 

c) Each column y, for which the elements of rows a and f are 
both greater than or equal to one, contains no finite negative numbers. 

d) The fractional elements are identically placed in each col- 
umn, and if there are arbitrarily large negatives, these occur in pairs, 
except that one column may contain additional negatives if it con- 
tains no elements greater than or equal to one. 

e) Neither column a nor column # contains arbitrarily large 
negative elements unless they occur in the same rows for both col- 


umns. 
f) Either column « or column / contains no elements which are 


proper fractions. 
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Then, if conditions (a), (b) and (c) hold, together with either 
condition (d) or conditions (e) and (f), one can delete one row and 
its corresponding column in F and the corresponding space in a if 
the two columns are combined in such a way that the algebraically 
larger of corresponding elements is retained in each case. 

It is very likely to be the case that the duplicated net of theorem 
(1) can be reduced to a smaller net which is still an expanded equiva- 
lent net of the original net. Suppose that the structure matrix of the 
original net contained no pair of rows or pair of columns which are 
equal. Then in the duplicated net, no two columns are equal, but there 
may well be rows which are equal. Suppose this is the case and a 
reduction is made. Then the corresponding F’; in the diagonal of 
F",, the corresponding section F’, and the corresponding section 
F’, of F” are altered. Suppose successive reductions are made only 
if pairs of identical rows pass through the same F”, in the diagonal 
of F”;. When there are no more obvious reductions of this kind left 
to be made, each of the non-zero sections of F'”, except F’”, , may have 
been altered. Thus we may write the following theorem (II): 

For any net of N’ neurons, N’ =p +1 + «, with a corresponding 
structure matrix F’, there is an extended equivalent net, a separated 
net, of N* neurons, N*=p +4 + ---+tat-:: te +e, which has @ 
structure matrix F* equal to F’x plus the extended sum of « matrices 
represented. by the (p + ta + €)? matrix F, which has zeros in all the 
last « columns except column p + ta + a, the first p elements of this 
column being zeros also. 


The resulting matrix F* is one in which F*, is given by diag 
{Fir, +++, Far, +++, Fer}, F*y = F’y, and F*; may be written as diag 
{6,,---,ba,---, bc}, where bc is a (ta X 1) column vector. 

The above theorem reduces the problem of the construction of a 
net with « efferents to the problem of the construction of a net with 
one efferent. Having constructed a net for each efferent alone, the 
total net can then be determined by theorem (II). We shall find, how- 
ever, a method of constructing nets which will not require use of the 
above theorem. We pass now to this problem: _ 

Suppose that a certain sequence Rao(r) of r’s leads to a finite 
sequence €o(e) of e’s, and that a number of such pairs of sequences 
is known empirically. It will be assumed that the net is always at 
rest when any afferent sequence is initiated. We wish to construct a 
net for which the structure matrix F(R, &) is such that each R, of 
a set ‘R of sequences of r’s leads to the appropriate sequence 6 of 
e’s. Suppose that ec(1) is the first non-zero e. If a particular se- 
quence R,(r) leads to a particular &.(e), we can write this sequence 
of r’s, ro(—l’), --- , ro(—r), «++ , Po (0), «++ , ro(l”) which results in a 
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sequence €c(1), €c(2), --- , ec(t),--- , es (l’”), as a set of sequences 
in which ro(—l’), ---,ro(0) results in es(1), and re(—l'), «++ , re(1) 
results in es(2), and so forth. Now let 8, be the first sequence, Ra 
be sequence o , and ®, be the last sequence, all sequences having been 
ordered. Let the time be adjusted for each sequence so that the last 
member of R is r(0), and each sequence &(e) contains but one ele 
ment e(1). Let i + 1 be the longest series of r’s so that, in general, 
sequence ‘Ro will begin with ro(—l) and end with r.(0), and have 
the general term ro(—7). 

Now ro(—r) will always be composed of zeros and ones, but some 
of these may not be known. Thus there will be some number mu-, 
ones, some number of zeros and a number of unknowns z. Similarly 
€c(1) will contain elements which may be one, zero or # . 

For the following construction, the occurrence of a one in the 
r-th element of rc(—7) implies that the 7-th efferent must act at time 


¢ = —; in order that e.(1) have its given value. The occurrence of 
a zero means that the corresponding afferent must not act at time 
t = —; if eo(1) is to have the value given. The occurrence of an x 


implies that it is immaterial whether or not the corresponding ele- 
ment acted at t = —-. It is understood that all r(¢) vectors for t < —l 
contain only x’s. In the case of the efferent vector, e , the occurrence 
of a unity in the e-th element implies that the sequence Ra does cause 
the e-th efferent to act at t = 1, whereas a zero in this element im- 
plies that sequence Ra prevents its action even though some part of 
R, could otherwise produce activity. But an z in eo(1) means that 
it is immaterial whether or not this element occurs as a result of R., 
though generally it may be considered that R. simply fails to produce 
activity in the corresponding efferent at?=1. 

Define an operator Bp which acts on a (1 X p) afferent activity 
vector ro(—7) such that Be{rco(—7)} is a (p X e) matrix made up of 
equal column vectors whose elements are obtained from ro(—7) in 
the following manner. Let n. be the number of vectors in Re which 
contain at least one unity. An element one in ro(—7) is changed to 
1/morNc , an element x is changed to zero, and an element zero is 
changed to —é where 6 is an arbitrary positive quantity. If rc (—7) 
contains at least one unity, then ro(—r) Br{ro(—7’)} is a row vector 
of repeated elements 1/nc if ro(—r) = ro(—7’), otherwise each ele- 
ment is less than 1/n.. When ro(—7) contains no ones, then the ele- 
ments of ro(—r) Be{ro(—7’) } are all zero if ro(—7) = Yo (—7'), other- 
wise each element is less than or equal.to zero. In a somewhat similar 
manner, define B; as an operator which acts on a (1 X e) efferent 
activity vector e such that By{e} is a (e X «) diagonal matrix whose 
diagonal elements are obtained from e by changing zero to —dé, «x to 
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zero, ones remaining unchanged. 

We shall next define an operator which extends a structure matrix 
by adding an element between each afferent and the neurons being 
acted upon by afferents, displacing these connections. If F is any 
structure matrix, then the afferent extending operator T, is defined 
by the relation: 


p 
ofa Pe etc 
If F=0 F, F,, then T;F=> R Fees (2) 
0 GeO 0. OF lz 
Ose Opes 0 


where 1° is a unit matrix of order p. The operator Tz is distributive 
with respect to extended sums if matrices which represent equivalent 
nets are said to be equal. 

Let F” be a structure matrix such that the sequence Ro less the 
last vector ro(0), which we assume to contain no ones, results in 
éo(0) = (1,1,1---, 1); thatis, e.(1) occurs one unit of time sooner 
than it should, eh anah that every efferent acts. Then it can be seen 
that 

Fg = X{ro(—1) } + TeX {ro (—2) } (4+) --- (4+) Te! X{re(—l) } 
(3) 
=D Tp X{ro(—r)}, 
where X{ro(—r)} is a structure matrix, (p + ¢«) X (p + «), which 
contains Br{ro(—r)} in its upper right hand corner but is elsewhere 
zero. The summation sign with a superscript (F£) will be used to 
denote an extended sum. 

Now F-. will contain a submatrix Fc, in which each column con- 
tains Be{ro(—1)}, Be{ro(—2)}, ---, and Br{ro(—l)}. Thus each ef- 
ferent will act as a result of a sequence ®’, if, and only if, the sum 
r’o(—1) Be{ro(—1)} + --- + 9r'o(—l) Be{ro(—l)} is greater than or 
equal to unity in any element. Now this sum cannot exceed %o(1/nc) 
and will equal unity if, and only if, r’s(—r) = ro(-<r) for every 7; 
hence, only the sequence Ra is effective. 

Analogous to the operator Tz, we can define an operator T,, 
which extends the structure on the efferent side, by 


(ren ee eo 
ie peat so: 4 
Ok Sah eae te (4) 
Teme cae 


the quantity 1° being the unit matrix of rank e. 


Ti Se 
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From the fourth rule for reducing a matrix it follows that 
Ta(F" (+) F’) = TF’ (+) T,F’ if Fx, F"y, F’;, F's contain no 
proper fractions or negative elements. Evidently, then, if F’ = T,F 
and F” = T,F,, Ty, will be distributive with respect to F’ and F’”. 

It will be necessary to use a more special operator T, which acts 
on a structure matrix as above but in which 1° is replaced by the 
diagonal matrix Bz{e}. Then T.F’s is a structure matrix such that 
the sequence Ra. , and only this sequence [r.(0) containing no ones], 
results in e(1), with the exception that ro(0) has no effect. Let 
X{ro(0), €c(1)} be a matrix of suitable order, which is zero every- 
where except in the last e columns of the first ¢ rows, this non-zero 
submatrix being [Br{ro (0) }] [Be{ec(1) }] . This matrix takes into ac- 
count the effect of the vector r.(0). Thus adding this to T, F’s results 
in a matrix F, representing a net in which the afferent sequence Ra 
results in the efferent activity vector es(1). The extended sum of 
these matrices represented by F. is then the structure matrix desired, 
and we have the following result: 

If each of a set R of sequences Ro of afferent activity actors re- 
suits in a corresponding sequence 6. of afferent activity vectors of a set 
€, and these have been written in such a way that each afferent se- 
quence contains only one vector, and if the last vector of each sequence 
of R contains no ones or if we choose to ignore them, and if e.(1) con- 
tains no zeros whenever ro(0) contains zeros, then the structure 
which converts a sequence of ® into a sequence & is represented by 
the matrix F(R, €) given by 


F(R, 6) = > X{r.(0), €c(1)} + Se) Te 


O=1 C=} 
I (5) 
wh) Ts." X {Tez} 315 
T=1 
the operators and symbols having been defined above. If no e con- 
tains a zero, then we may also write 


F(R, 6) = DX {ro(0) , €o(1)} + =o Ts 


3 O=1 O=1 (6) 
AD Ons X{Yed—7),, €0(1) 31; 


br igs I 


If the last activity vectors contain ones, then the above method 
can be modified to give a result under certain conditions. If it is the 
ease that in forming the sum of X{ro(0), ec(1)} over o, no element 
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is the result of adding a positive term with either a negative or a posi- 
tive term, unless both are greater than or equal to one, and if no €o(1) 
contains zeros, then F(R, €) can be written as in equation (6) but 
with Tz replaced by Tp. Since Tr(Fi (+) Fs) = TeFi (+) TeF after 
reducing by rule (2), the operator Tz outside the brackets can be 
brought under the second summation sign. Since an ordinary sum in 
the first term of equation (6) cannot be confused with the extended 
sum, we may then write for this case 


F(R, €) =S" S Ty X(re(—2) , eo (1)}. (7) 


° T=0 


More generally, however, when ro(0) contains ones in a certain 
subgroup of sequences, o=s’' + 1,---,s, then instead of equation (7) 
we may write an equation in which the sum over o from 1 to s’ is giv- 
en as in equation (5), while the sum over o from s’ + 1 to s is given 
as in equation (7). This gives a more general equation for the con- 
struction of a net: 


F(R, 6) =>. X {ro(0) , €0(1)} 


O=1 


(+) SH oe x Te AAT ny (8) 


O=1 T=1 


8 1 
(+) DD Ta X (re(—z) , eo(1)}. 
O=8'+1 T=0 

If no e(1) contains a zero, T, may be replaced by Tz in equation (8): 

Let the superposition of two sequences be a sequence formed by 
combining corresponding elements (1, x or 0) of corresponding vec- 
tors as follows: The combination of an element with itself is the ele- 
ment itself. The combination of an x with any element is that element © 
with which x is combined. If the sequence is an efferent, then the com- 
bination of a zero and a one is the element zero. If the sequence is an 
afferent, then if zero and one combine in any element, the superposi- 
tion is a null sequence. The superposition of sequences is commuta- 
tive and associative. 

Let a set of afferent and efferent sequences (R, €) be given. Let 
(R, €)r be the set containing (R, €) and all its non-null superposi- 
tions. Then we may refer to a net as being structurally equivalent to 
a set of sequence pairs (R, €) if every Rein (R,&), applied to the 
net results in &¢ but every other afferent sequence results in 6 = 0. 
If the above holds for every pair in (R, €) but not for all superposi- 
tions, then we may refer to the net as being structurally equivalent to 
(R, €) ina stronger sense. Similarly, we may refer to a net as being 
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adequate with respect to (R, €) if each given afferent sequence ap- 
plied to the net results in its corresponding efferent sequence but some 
other sequences, not in ® , result in non-null efferent sequences. 

If equation (8) reduces to equation (5) or (6), then the net con- 
structed is structurally equivalent to (R, €). The net obtained from 
equation (8) may be structurally equivalent even though that con- 
structed by equation (7) is only adequate; but the net obtained from 
equation (7) may be equivalent in a stronger sense. The net obtained 
from equation (8) will be equivalent, perhaps in the stronger sense, 
if the sum of the elements of every efferent column of the structure 
matrix which contains proper fractions is unity. 

Under the following conditions it will not be possible by the above 
method to construct a net which is even structurally adequate. Let 
o,o,7, and e be subscripts which can take on any value in their re- 
spective ranges [(1, ---,s),(1,---, s), (1,>--, p), (1,-:-, 2)]. 
The conditions are then: a) ro(0) and e.(1) both contain zero ele- 
ments, and b) the scalar elements 7,(0) = €ce(1) = 1, while either 
fo,(0) = 1 and ec. (1) = 0 or 75',(0)=—0 and eg (1) =1. 

In constructing a net, one may use the rules for reducing nets 
but not the rules for changing elements. Since the quantity —6é is of 
arbitrary magnitude, one may choose it to be arbitrarily large. In 
matrices, this will be indicated simply by a minus sign. By choosing 
6’s large, the probability of being able to reduce the net is increased. 
It should also be noted that the use of the operator T, in equation (5) 
results in ¢ identical columns of which all but one can be eliminated 
if the corresponding columns are added. This results in a row vector 
in place of the diagonal matrix, the elements of the row being the 
same as the elements of the original diagonal submatrix. 

The construction of a net from equation (8) may be illustrated 
by an example from W. S. McCulloch and W. Pitts (1943) cited pre- 
viously (Landahl and Runge, 1946). In this case we have R, = [(2,1), 
(x, 0); (wx)J,e. = (1,2); R2= (1, %), & = (1, #) and R; = 
[(z, 1), (w, 1)], e: = (@, 1). Readjusting the time so that each e 


is e(1), then Ba(r.(0)}=0, Be(ri(—1)} = [2°], Belrs(—2)) = 


00 
LT | Belo) = ool > Bers} = [4 | > Bel) 
= Lo y Similarly Br {ei} = Br{e2} es ; ; and Br {es} 
j2 2 : 
00 : = 1 0 
= ir 1h" The Fy section of the sum X{ra (0), @o(1) } =0 + 00 


i eae oy! 
|0 4 04 
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Now Tr X {r,(—2)} is a (5 X 5) matrix whose F; , F, and F,; 
sections are, respectively, : , {1 1| and zero. Thus the term due to 


sequence R,, Te [X{ri(—1)} (+) Te X{ri:(—2)}], is a (7 X 7) ma- 
trix which reduces to the (6 X 6) matrix for which the sections Fz , 


| 1 ‘ 
F, and F; are, pas meeeN NE : - : | : ; ; 0 . The last term 


Tz X {rz(—1), e3(1)} is a (5 X 5) matrix whose Fz, Fr and F, sec- 


, 


tions are ? | , |0 4|, and zero. The extended sum of the above 
matrices leads to a (7 X 7) matrix in which there are two identical 
columns. Elimination of one row and column leads to a matrix which 
is identical to that previously written (Landahl and Runge, 1946) to 
represent the net of W. S. McCulloch and W. Pitts (1943), except for 
the negative element. In this case this element is arbitrary, while by 
the other construction it is an arbitrarily large negative, though it 
was written as minus one for convenience. It is evident by inspection 
that the magnitude of this element does not affect the equivalence of 
the nets in this case. The-resulting net is structurally equivalent to 
the given set of pairs of sequences. 

If, instead of using equation (8): for constructing the net, equa- 
tion (7) had been used, then the resulting matrix would differ from 
that obtained above in that the element of row 2 and column 4 would 
be zero instead of negative, while the element of row 3 column 5 
would be negative instead of zero. This net is structurally equivalent 
in the stronger sense to the given set of pairs of sequences. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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